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ABSTRACT 

We present two independent, homogeneous, global analyses of the transit lightcurves, radial veloc- 
ities and spectroscopy of Kepler-4b, Kepler-5b, Kepler-6b, Kepler- 7b and Kepler-8b, with numerous 
differences over the previous methods. These include: i) improved decorrelated parameter fitting set 
used, ii) new limb darkening coefBcients, iii) time stamps modified to BJD for consistency with RV 
data, iv) two different methods for compensating for the long integration time of Kepler LC data, v) 
best-fit secondary eclipse depths and excluded upper limits, vi) fitted mid-transit times, durations, 
depths and baseline fluxes for individual transits. We make several determinations not found in the 
discovery papers: i) We detect a secondary eclipse for Kepler- 7b of depth (47 ± 14) ppm and statistical 
significance 3.5-it. We conclude reflected light is a much more plausible origin than thermal emission 
and determine a geometric albedo of Ag = (0.38 ± 0.12). ii) We find that an eccentric orbit model for 
the Neptune-mass planet Kepler-4b is detected at the 2-a level with e = (0.25 ± 0.12). If confirmed, 
this would place Kepler-4b in a similar category as GJ 436b and HAT-P-llb as an eccentric, Neptune- 
mass planet, iii) We find weak evidence for a secondary eclipse in Kepler-5b of 2-cr significance and 
depth (26 ±17) ppm. The most plausible explanation is refiected light caused by a planet of geometric 
albedo Ag = (0.15 ± 0.10). iv) A 2.6-cr peak in Kepler-6b TTV periodogram is detected and is not 
easily explained as an aliased frequency. We find that mean-motion resonant perturbers, non-resonant 
perturbers and a companion extrasolar moon all provide inadequate explanations for this signal and 
the most likely source is stellar rotation, v) We find different impact parameters relative to the discov- 
ery papers in most cases, but generally self-consistent when compared to the two methods employed 
here, vi) We constrain the presence of mean motion resonant planets for all five planets through 
an analysis of the mid-transit times, vii) We constrain the presence of extrasolar moons for all five 
planets, viii) We constrain the presence of Trojans for all five planets. 

Subject headings: planetary systems — stars: individual (Kepler-4, Kepler-5, Kepler-6, Kepler-7, 
Kepler-8) techniques: spectroscopic, photometric 



1. INTRODUCTION 

The Kepler Mission was successfully launched on 
March 7th 2009 and began science operations on May 
12th of the same year. Designed to detect Earth- 
like transits around Sun-like stars, the required pho- 
tometric precision is at the level of 20 ppm over 6.5 
hours integration on 12"^ magnitude stars and early re- 
sults indicate this im pressive precision is being reached 
(jBorucki et all I2009D . In 2010, the first five transit- 
ing exoplanets (TEPs) to be discovered by the Kepler 
Mission were announ ced by the Kepler Science Team 
(|Borucki et al.l I2010a[) . known as Kepl er-4b, Kepler- 
5b, Kep le r-6b, Kepler - 7 b and K epler-8b (iBorucki et al.l 
(I2010bl): iKochct al.' (20To); Dunham^^alJ ([20100; 
iLatham et al. (2010); Jenkins et al. (2010at)). These ex- 
panded the sample of known transiting exoplanets to 
about 75 at the time of announcement. 

The main objective of the Kepler Mission is to dis- 
cover Earth-like planets, but the instrument naturally 
offers a vast array of other science opportunities includ- 
ing detection o f gas giants, searches for thermal emission 
(|Deming et al.li2005.) and/or reflec tion from exoplanets , 
detection of orbital phas e curves (iKnutson et al.ll2007| ) 



and ellip soidal variation s ('Wels h et ah 
mology (|Christensen-Da lsgaard ~t al.l 
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timing (|Agol et all [20051 ). to name a few. Confirma- 
tion and follow-up of exoplanet transits is known to be 
a resource intensive activity and since the detection of 
new planets is Kepler's primary objective, it is logical 
for many of these other scientific tasks to be conducted 
by the astronomical community as a whole. 

Independent and detailed investigations of the Kepler 
photometry provides an "acid-test" of the methods em- 
ployed by the Kepler Science Team. Indeed, the distinct 
analysis of any scientific measurement has always been a 
fundamental corner stone of the scientific method. In this 
paper, we present two independent analyses of the dis- 
covery photometry for the first five Kepler planets. We 
aim to not only test the accuracy of the methods used in 
the discovery papers, but also test our own methods by 
performing two separate studies. Both methods will be 
using the same original data, as published in the discov- 
ery papers. 

Some additional data-processing tasks are run through 
the Kepler reduced data, which we were not used in the 
original analyses presented in the discovery papers, and 
are discussed in §3. In this section, we also discuss the 
generation of new limb darkening coefficients and meth- 
ods for compensating for the long integration time of the 
Kepler long-cadence photometry. We also perform indi- 
vidual transit fits for all available Kepler transits in or- 
der to search for transit timing variations (TTV), transit 
duration variations (TDV) and other possible changes. 
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These will be used to provide a search for perturbing 
planets and companion exomoons. 

2. FITTING METHODOLOGY 

2.1. Method A 

In this work, the two of us adopt different algorithms 
for fitting the observational data. The first code has been 
written by D. Kipping and is a Markov Chain Monte 
Carlo (MCMC) code (for a bri ef introduction on th is 
method, consult appendix A of iTegmark et al.l (|2004[ )) 
with a Metropolis-Hastings algorithm, written in For- 
tran 77/90. The lightcurve mode l is generated using the 
analytic quadratic limb darkening lMandel &: AgoU ()2002[ ) 
routine, treating the specific stellar intensity with: 



= 1-ui(1-m)-M2(1-m)' 



(1) 



The quadratic limb darkening coefficients_wi and U2 
are known to be highly correlated (jPall I2008bl) . but a 
principal component analysis provides a more efficient 
fitting set: 



u'l = ui cos 40° — U2 sin 40° 
1*2 — ui sin 40° 4- U2 cos 40° 



(2) 
(3) 

We also implement conditions to ensure the brightness 
profile is everywhere positive and monotonically decreas- 
ing from limb to c enter: ui > and < ui -I- U2 < 1 
(jCarter et al.l [20?)8h. 

Quantities relating to time differences are computed 
by solving the bi-quartic equation for the various true 
anomalies of interest followin g the m ethods detailed in 
IKippingj (|2008l ) and IKippini3 (|2010aD . These quantities 
include: 

• The transit duration between the first and fourth 
contact, Ti,4. 

• The transit duration between the second and third 
contact, T2,3. 

• The ingress and egress transit durations, T12 and 
T34 respectively. 

• The secondary eclipse full duration, S'1,4. 

• The time of the predicted secondary eclipse, igcc- 

• The offset time between when the radial velocity 
possesses maximum gradient (for the instance clos- 
est to the time of conjunction) and the primary 
mid-transit, AiRy. 

The latter two on this list are particularly important. 
For example, the RV offset time provides a strong con- 
straint on ecosw. 

The principal fitting parameters for the transit model 
are the orbital perio d, P, the ratio-of-radii squared, p^, 
the IKippind (|2010af) transit duration equation, Ti , the 
impact parameter, b, the epoch of mid-transit, E and 
the Lagrange eccentricity parameters k = e cos w and 
h = esinw. Ti is used as it describes the duration be- 
tween the planet's centre crossing the stellar limb and 
exiting under the same condition and this is known to 



be a useful decorr elated parameter for light curve fits 
(jCarter et al.l 120081 ). We also note that an approximate 
form is required as the exact expression requires solv- 
ing a bi-quartic equation which is not easily invertible. 
In contrast, Ti may be inverted to provide a/i?* using 
(|Kippindl20T(M ): 



Qc = 



e sm uj 



(4) 



a/R., 



\ 




Blendin g is acco unted for using the prescription of 
iKipping fc Tinettil (2010) where we modify the formu- 
las for an independent blend rather than a self-blend. 
The nightside pollution effect of the planet is at least 
an order of magnitude below the detection sensitivity of 
Kepler due to the visible wavelength bandpass of the in- 
strument. Therefore, the only blending we need account 
for are stellar blends. In addition, the model allows for 
an out-of-transit fiux level (Foot) to be fitted for. Unlike 
method B (see later), we fix Foot to be a constant during 
the entire orbital phase of the planet. S econdary eclipses 
are produced using the lMandel fc Agoll (|2002) code with 
no limb darkening and the application of a transforma- 
tion onto the secondary eclipse lightcurve which effec- 
tively squashes the depth by a ratio such that the new 
depth is equal to {Fpd/F^) (planet dayside flux to stel- 
lar fiux ratio) in the selected bandpass. This technique 
ensures the secondary eclipse duration and shape are cor- 
rectly calculated. The observed fiux is therefore modeled 
as: 



F, 



obs 



fmodel(t) + (^ - 1) 

B 



Fo, 



(6) 



This model so far accounts for primary and secondary 
eclipses but an additional subroutine has been written 
to model the radial velocity variations of a single planet. 
Modeling the RV in conjunction with the transit data 
gives rise to several potential complications. Firstly, we 
could try fits using either an eccentric orbit or a fixed 
circular orbit model. An eccentric fit will always find a 
non-zero eccentricity du e to t he boundary condition that 
e > (Lucy & Sweeney 1971f ) . The effects of fitting ver- 
sus not-fitting for eccentricity are explored in this work 
by presenting both fits for comparison. 

A second complication is the possible presence of a 
linear drift in the RVs due a distant massive planet. 
This drifts can give rise to artificial eccentricity if not 
accounted for but their inclusion naturally increases the 
uncertainty on all parameters coming from the RV fits. 
Thirdly, an offset time, At^v, between when then RV 
signal has maximum gradient (for the instance nearest 
the time of conjunction) and the mid-transit time can 
be included. For eccentric orbits, a non-zero offset al- 
ways exists due to the orbital eccentricity. In our code 
this is calculated and denoted as Atecc and is calculated 
exactly by solving for the relevant true anomalies and 
co mputing t he tim e interval using the duration function 
of IKipping ()2008[ ). However, an additional offset can 
also exist, Aftroji which may due to a massive body in a 
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Trojan orbit (i.e. we have A^rv = A^qcc + A tny). This 
offset was first predicted by iFord fc Holmanl (2007) and 
the detection of a non-zero value of Aitroj would indicate 
a Trojan. For eccentric orbits, the error in A^occ is of- 
ten very large and dominates the error budget in A^troj j 
washing out any hints of a Trojan. 

The question exists as to when one should include 
these two additional parameters. Their perenial inclu- 
sion would result in very large errors for poorly charac- 
terized orbits and this is clearly not desirable. We there- 
fore choose to run a MCMC -I- simplex fit for all possible 
models and extract the lowest ^ for the RV signal. This 
minimum is then used to co mpute th e Baye s ian In- 
forma tion Criterion, or BIG fsee iSchwarzl ()1978[ ): iLiddld 
(|2007f) ). BIG compares models with differing numbers of 
free parameters (fc), heavily penalizing those with more 
and the preferred model is given that yielding the lowest 
BIG, where: 

BIG = + fclogA^ (7) 

Where iV is the number of data points. In total, there 
exists four possible models for the circular and four pos- 
sible models for the eccentric fit by switching on/off these 
two parameters. For both the circular and eccentric fits, 
we select the model giving the lowest BIG. Therefore, if 
for example the lowest BIG was that of a zero 7 but non- 
zero Attrojj we would fix the gradient term but let the 
offset be freely fitted in the final results. The results of 
these preliminary investigations can be found in Table [T^ 
of the Appendix. 

The favored solution is generally the eccentric model 
over the circular model. This is because the light curve 
derived stellar density, which we will later use in our 
stellar evolution models, has a strong dependence on the 
eccentricity. In general, fixing e = leads to unrealisti- 
cally small error on p*. However, in some cases sparse 
RV phase coverage can lead to artificially large e values. 

The code is executed as a global routine, fitting all the 
RV and photometry simultaneously with a total of > 13 
free parameters: -E, p^, Ti, 6, P, itj, h — esincj, 
k = ecosw, 7rci, Foot and (Fpd/F-i,). The inclusion 
of 7 and A^troj is reviewed on a case-by-case basis, as 
discussed. Additionally, the blending factor, which we 
denote as B, is allowed to float around its best fit value 
in a Gaussian distribution of standard deviation equal to 
the uncertainty in B. These values are reported in the 
original discovery papers. 

The final fit is then executed with 125,000 trials with 
the first 20% of trials discarded to allow for burn-in. This 
leaves us with 100,000 trials to produce each a poste- 
riori distribution , which is more than adequate. The 
IGelman &: RubinI ([1992) statistic is calculated for all fit- 
ted parameters to ensure good mixing. The final quoted 
values are given by the median of the resultant MGMG 
trials. We define the posian and negian as the maxi- 
mum and minimum confidence limits of the median re- 
spectively. The negian may be found by sorting the list 
and extracting the entry which is 34.13% of the total list 
length. The posian is given by the entry which occurs 
at 68.27% of the total list length. These values may be 
then be used to determine the confidence bounds on the 
median. 

After the global fitting is complete, we produce the dis- 



tribution of the lightcurve derived stellar density and a 
Gaussian distribution for the SME-derived effective tem- 
perature and metallicity around the values published in 
the discovery papers and based on SME analysis of high 
resolution spectra. These distributions consist of 100,000 
values and thus may be used to derive 100,000 estimates 
of t he stellar prop erties through a YY-isochrone analy- 
sis (lYi et al.l[200l . For errors in TcSi, and [Fe/H], we 
frequently adopted double that which was quoted in the 
discovery papers, as experience with HAT planets has re- 
vealed that SME often underestimates these uncertain- 
ties. Finally, the planetary parameters and their uncer- 
tainties were derived by the direct combination of the a 
posteriori distributions of the lightcurve, RV and stellar 
parameters. The stellar jitter squared is found by tak- 
ing the best-fit RV model residuals and calculating the 
variance and subtracting the sum of the measurement 
uncertainties squared. 

2.2. Method B 

The second method builds on the codes developed un- 
der the HATNet project, mostly written in G and shell 
scripts. These have been used for th e analys i s of HATNet 
planet disco veries, s uch as lPal et"aLl (j2008al) : rBakos et all 
(j2009[ ):|Plil (j2009bD . and have been heavily modified for 
the current case of Kepler analysis. We used a model for 
describing the F{t) fiux of the star plus planet system 
through the full orbit of the planet, from primary transit 
to occultation of the planet. In our initial analysis we as- 
sumed a periodic model, and later we considered the case 
of variable parameters as a function of transit number. 
The F{t) fiux was a combination of three terms: 

F{t) = E, P, Foot, Foos,Foo,n, C/R*,P, b\ k, h, w) - 

Foot/{l + epn + B)x Ft{t, E, P, C/i?*,P, &') - 

Foos X €pd/{{l + €pd + B)p^) 

xFs{t,E,PX/R.,P,b'') (8) 

Here $(.) (the first term) is the brightness (phase- 
curve) of the star+planet-f blend system without the ef- 
fect of the transit and occultation. The drop in flux due 
to the transit is reflected by the second term, and the 
equivalent drop due to the occultation of the planet by 
the third term. The $(.) phase-curve is a rather sim- 
ple step-function of three levels; Foot within w times the 
duration of the transit around the transit center, Foos 
within w times the duration of the occultation around 
the center of the occultation, and Foom in between. For 
simplicity, we adopted w = 2. Because the transit and 
occultation centers and durations also enter the formula, 
$0 depends on a number of other parameters, namely 
the t time, E mid-transit time, P period, the C/R* Pa- 
rameter that is characteristic of the duration of the tran- 
sit, p = Rp/R-i, planet to star radius ratio, impact 
parameter, and k and h Lagrangian orbital parameters. 
The location and total duration of the transit and occul- 
tation are fully determined by the above parameters. In 
general, the combined brightness Foom of the star plus 
planet would smoothly vary from Foot outside of the pri- 
mary transit to Foos outside of the occultation, typically 
as a gradual brightening as the planet shows its star-lit 
face towards the obse rver. If the night-side brightness of 
the planet is ep„ • F^, ([Kipping &: Tincttl.2010,) . and the 
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day-side brightness of the planet is epd ■ i^*, and BF^, is 
the flux of a blend contributing light to the system, then 
1 + B + Cpn = Foot and 1 + B + epd = Foos (both equations 
normalized by Fi,), i.e. there is a constraint between Foot 
and Foos- The reason for introducing all three of footj 
Foom and Foos for characterizing the data was because 
the Kepler lightcurve, as provided by the archive, was 
"de-trended" , removing the (possible) brightening of the 
phase-curve due to the planet. Because of this, we found 
that introducing a more complex phase function was not 
warranted. We note that while such de-trending may be 
necessary to eliminate systematic effects, it also removes 
real physical variations of small amplitude. 
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Fig. 1. — Top: The phase-curve of our simplistic model, show- 
ing a primary transit, an occultation, and a fix flux level Foom 
in between. The model was generated with unrealistic planetary 
parameters to aid the visualization: Rp/ Fd, = 0.2 (large planet-to- 
star radius ratio), = o.25, C/R* = 10, B = 0, P = 3.0, k = 0.1, 
h = -0.2 (eccentric orbit), ep„ = 0.01, tp^j = 0.02, B = 0.2, 
Foot = 1.2, Foos = 1.3, Foom = 1.25 and w = 2.0. See text for the 
meaning of these parameters. Due to the eccentric orbit, the occul- 
tation is clearly off'set from halfway between the primary transits, 
and due to the lack of limb-darkening, the occultation exhibits a 
flat bottom. 

The second term characterizes the dimming of the com- 
bined light of the star, blend and planet during the pri- 
mary transit. The loss of stellar flux Ft(.) was mod- 
eled u sing the analytic formula based on lMandel &: Ago! 
(I2002D for the eclipse of a star by a planet (|Bakos et aL 
120091) ■ where the stellar flux was described by quadratic 
limb-darkening. The quadratic limb-darkening coeffi- 
cients were derived by evaluating the stellar flux as a 
function of = cos 6 (where 6 is the angle between the 
stellar surface normal vector and the line-of-sight) us- 
ing the quadratic ui, U2 coefficients. The formula takes 
into account the B blending (constant flux contribution 
by other sources) , the ep„ night-side illumination of the 
planet (normalized by F^,), and the Foot out-of-transit 
flux level that was normalized by an arbitrary constant 

Finally, the third term characterizes the dimming of 
the combined light of the star, blend and planet during 
the occultation of the planet. The day-side flux of the 
planet Fpd is gradually lost as it moves behind the star. 
This is modeled by a zero limb-darken ing model Fs(.) 
of the secondary eclipse using the same IMandel fc Asol 
(|2002l l formalism. This term also takes into account the 
B blending, and the Cpd daytime flux of the system. 

^ Had the flux been normalized by the Fi, brightness of the star, 
without blending and other sources. Foot = 1 would hold. 



Our global modeling also included the radial velocity 
data. Following the formalism presented by iPall ([2009) , 
the RV curve was parametrized by an eccentric Keple- 
rian orbit with semi-amplitude K, Lagrangian orbital el- 
ements {k,h) = ex (coso;, sincj), and systemic velocity 

7rcl- 

The Kepler photometry and radial velocity data are 
connected in a number of ways. For example, we as- 
sumed that there is a strict periodicity in the individual 
transit time^, and the same E and P ephemeris of the 
system describes the photometric and RV variations. An- 
other example: the Lagrangian orbital parameters k and 
h were not only determined by the RV data, but also by 
the phase of the occultation of the planets in the Kepler 
data. 

Altogether, the 13 parameters describing the physical 
model were Ta (the Tc of the first transit in the Ke- 
pler data), Tb (the same parameter for the last tran- 
sit), K, k = ecosw, h = esinw, 7rei, Rp/R^,, 6^, 
C/i?*, epd, Foot, Foom, and Foos- The B blending fac- 
tor was kept fixed a t the values published in the Ke- 
ier di sc overy papers (iBorucki et al.l (l2010bD:lKoch et al.l 
20T0l): iDunham et al.l (|2010D : iLatham et al.l (|2010D : 



iJenkins et al.l ()2010aD ). We adopted ep„ = for the 
night-side emission of the planet. Fitting for B and ep„ 
would require data with exceptional quality. 

The joint fit on the photometry and radial velocity 
data was performed as described in [Bakos ct al. (2009). 
We minimized in the parameter space by using a hy- 
brid algorithm, combining the do wnhill simplex method 
(AMOEBA; see lPress et"allll992[ ) with the classical lin- 
ear least squares algorithm. Uncertainties for the param- 
eters were derived usi ng the Mark ov Chain Monte-Carlo 
method (MCMC, see [S3 [IMl). The a priori distri- 
butions of the parameters for these chains were chosen 
from a generic Gaussian distribution, with eigenvalues 
and eigenvectors derived from the Fisher covariance ma- 
trix for the best-fit value. The Fisher covariance matrix is 
calculated analytically using the partial derivatives given 
by Pal (200j). 

Stellar parameters were also determined in a Montc- 
Carlo fashion, using the t ools and methodology described 
in e.g. lBakos et~all ([20091) . Basically, we used the MCMC 
distribution of a/Ri, and corresponding p^,, a Gaussian 
distribution of Tcff^, and [Fe/H] around the values pub- 
lished in the discovery papers and based on SME analysis 
of high resolution spectra. For errors in Teg* and [Fe/H] 
we typically adopted double that of quoted in the discov- 
ery papers. Finally, the planetary parameters and their 
uncertainties were derived by the direct combination of 
the a posteriori distributions of the lightcurve, RV and 
stellar parameters. 

Throughout this text, the final value for any given 
parameter is the median of the distribution derived by 
the Monte-Carlo Markov Chain analysis. Error is the 
standard deviation around the final value. Asymmetric 
error-bars are given if the negative and positive standard 
deviations differ by 30%. Parameter tables also list the 
RV "jitter," which is a component of assumed astrophys- 
ical noise intrinsic to the star that we add in quadrature 

Note that we also performed a global analysis under method B 
by allowing the parameters of the individual transits (such as the 
transit centers) to vary. 



Kepler-4b through Kepler-8b 



5 



to the RV measurement uncertainties in order to have 
/{n — dof ) = 1 from the RV data for the global fit. 

2.3. Parameter variation search method 

2.3.1. Individual fits 

In addition to global fits, we will search for variations 
in the lightcurve parameters for each transit. Searching 
for changes in the mid-transit time, commonly known as 
transit timing variat ions (TTV) , can be u sed to search for 



pertu rbing planets (lAgol et al. 12005: H olman fc Murray! 
2005^ or companion moons (jSartoretti fc SchneiderillQQa 
Kipping 2009a). Similarly, transit duration variations 
(TDV) provid e a complemen tary method of searching 
for exomoons (|Kippindl2009allH ). We also look for depth 
changes and baseline variations which may aid in the 
analysis of any putative signals. 

Transit timing and duration measurements are made 
by splitting the time series into individual transits and 
fitting independently. In these individual fits, we only 
fit for , b, Ti, tc and Foot and float all other param- 
eters around their global best-fit determined value with 
their a posteriori distribution. The limb darkening coef- 
ficients are fixed to the theoretical values for these fits, 
since epoch to epoch variation in the LD coefficients is 
not expected. We choose to use method A for the in- 
dividual fits as the parameter b is more likely to have a 
uniform prior than b^ (as used in method B), from geo- 
metric considerations. As for the global fits, we perform 
125,000 trials for each fit with the first 25,000 (20%) tri- 
als being used as a burn-in. 

For the TDV, we define our transit duration as the 
time for the sky-projection of the planet's center to move 
between overlapping the stellar disc and exiting under 
th e same condi t ion, T . This value was first proposed 
by ICarter et al.l (j2008D due to its inherently low correla- 
tions with other parameters. As a result, this definition 
of the duration may be determined to a higher precision 
than the other durations. T has the property that its 
uncertainty is exactly twice that of the mid-transit time 
for a trapezoid approximated lightcurve. In reality, limb 
darkening will cause this error to be slightly larger. Nev- 
ertheless, this property means that a useful definition of 
TDV is to halve the actual variations to give TDV = 

0. 5(r - (T)) where (T) is the globally-fitted value of T. 
The factor of 0.5 means that we expect the r.m.s. scatter 
of the TDV to be equal to, or slightly larger than, the 
TTV scatter. It is also worth noting that in reality we 
could use any factor we choose rather than 0.5, because 
TDV is really a measurement of the fractional variation 
in duration and is not absolute in the sense that TTV is. 

2.3.2. Periodograms 

When analyzing the TTV and TDV signals, we first 
compute the of a model defined by a static system, 

1. e. constant duration and linear ephemeris. This 
value is naturally computed from the MCMC uncertain- 
ties for each transit. In practice, these errors tend to be 
overestimates. This is because the MCMC only produces 
accurate errors if it is moving around the true global min- 
imum. In reality, slight errors in the limb darkening laws, 
the finite resolution of our integration-time compensa- 
tion techniques and the error in the assumed photomet- 
ric weights themselves mean that our favorite solution 



is actually slightly offset from the true solution. Con- 
sequently, the MCMC jumps sample a larger volume of 
parameter space than if they were circling the true solu- 
tion. 

Despite the wide-spread use and therefore advocation 
of the MCMC method for transit analysis, it is also un- 
clear whether the MCMC method does produce com- 
pletely robust errors, particularly in light of the in- 
evitable hidden systemic effects. Therefore searching for 
excess variance in the values may not be a completely 
rel iable method of sea rching for signals, echoed recently 
bv lGibson et al.l (|2009( ). Concordantly, in all of our anal- 
yses, we compare the of a static model versus that 
of a signal and compute the F-test. The F-test has the 
advantage that it does care not about the absolute 
values, only the relative changes. Furthermore, it penal- 
izes models which are overly complex (Occam's razor). 

We therefore compute periodograms by moving 
through a range of periods in small steps of 1/1000 of 
the transit period and then fitting for a sinusoidal sig- 
nal's amplitude and phase in each case. We compute 
the x^ of the new fit and then use an F-test to com- 
pute the false alarm probability (FAP) to give the F-test 
periodogram. Any peaks below 5% FAP are investigated 
further and we define a formal detection limit of 3-cr con- 
fidence. Compared to the generalized Lomb-Scargle peri- 
odogram, the F-test is more conservative as it penalizes 
models for being overly-complex and therefore offers a 
more prudent evaluation of the significance of any sig- 
nal. Conceptually, we can see that the difference is that 
a Lomb-Scargle periodogram computes the probability 
of obtaining a periodic signal by chance whereas the F- 
test periodogram computes the statistical significance of 
preferring the hypothesis of a periodic signal over a null 
hypothesis. 

2.3.3. Analysis of variance 

In general, there are two types of periodic signals we 
are interested in searching for; those of period greater 
than the sampling rate (long-period signals) and those 
of periods shorter than the sampling rate (short-period 
signals). An example of a long-period signal would be 
an exterior perturbing planet and a short-period signal 
example would be a companion exomoon. Periodograms 
are well suited for detecting long period signals above the 
Nyquist frequency, which is equal to twice the sampling 
period. Short-period signals cannot be reliably seeked 
below the Nyquist frequency. First of all, a period equal 
to the sampling rate will always provide a strong peak. 
Periods of an integer fraction of this frequency will also 
present strong significances (e.g. 1/2, 1/3, 1/4, etc). As 
the interval between these aliases becomes smaller, very 
short-period signals become entwined. If a genuine short- 
period signal was present, it would exhibit at least partial 
frequency mixing with one of these peaks to create a 
plethora of short-period peaks. As a result, peaks in 
the periodogram below PNyquist cannot be considered as 
evidence for authentic signals. Therefore, we only plot 
our periodograms in the range of -PNyquist to twice the 
total temporal coverage. 

To resolve this issue, it is possible to derive an exo- 
moon's period by taking the rati o of the TTV to TDV 
r.m.s. amplitudes ('Kipp ins]|2009al lb[) . Since both of these 
signals are very short period, the data would ostensibly 
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be seen as excess scatter. These excess scatters can be 
used to infer the r.m.s. amphtude of each signal and then 
derive the exomoon's period from the ratio. Excess scat- 
ter, or variance, may be searched for by applying a simple 
test of the data, assuming a null-hypothesis of no sig- 
nal being present. If a signal is found, or we wish to 
place constraints on the presence of putative signal am- 
plitudes, we use the x^-distribution due to both Gaussian 
noise plus an embedded high frequency sinusoidal signal 
to model the variation. 

2.3.4. Phasing 

An unusual aspect of the individual data sets is that 
consecutive data points are 30 minutes apart. Despite 
accounting for the long integration times in our modeling 
(see later in §3.5), we consider the possibility that the 
long cadence data could produce artificial signals in the 
parameter variation search. We propose that one method 
of monitoring the effects of the long integration time is 
what we denote as phasing. 

Using the globally fitted ephemeris, we know the mid- 
transit time for all transits (assuming a non- variable sys- 
tem) . We may compare the difference between this mid- 
transit time and the closest data point in a temporal 
sense. This time difference can take a maximum value 
of ±15 minutes (assuming no data gaps) and will vary 
from transit to transit. We label this time difference as 
the phasing. In all searches for parameter variations, a 
figure showing the phasing will also be presented. This 
allows us to check whether any putative signals are cor- 
related to the phasing, which would indicate a strong 
probability the signal is spurious. 

3. DATA PROCESSING 

3.1. Barycentric Julian Dates 

The radial velocity data from Kepler is provided in 
barycentric Julian date (BJDutc) but the photometry is 
given in heliocentric Julian data (HJDutc)- For consis- 
tency, we calculate the difference between HJD and BJD 
for every time stamp and apply the appropriate correc- 
tion. We use the JPL Horizons ephemeris to perform 
this correction, which is important given that the differ- 
ence between HJD and BJD can be a few seconds. We 
do not correct times f rom the UTC system to TDB/TT 
(jEastman et al.l l201Cl[ ) as the difference is effectively a 
constant systematic due to leap seconds which accumu- 
late over decades timescale. Further, recent Kepler data 
products (e.g. Kepler Data Release 5) use the BJDutc 
system and thus there is a preference to remain consis- 
tent with the probable future data format. 

3.2. Measurement uncertainties 

Unfortunately, the Kepler pipeline output did not in- 
clude photometric uncertainties which forces us to eval- 
uate the measurement errors ourselves. For Kepler-4b, 
Kepler-5b, Kepler-6b and Kepler-7b, we found that the 
r.m.s. scatter was constant over the observational pe- 
riod. In order to estimate the standard deviati on, we 
use the median absolute deviation (MAD) from iGaussI 
11816) due to its robustness against outliers. Assum- 
ing Gaussian errors, the standard deviation is given by 
1.4286 multiplied by the MAD value. This value for the 
standard deviation is then assigned to all data points as 
their measurement uncertainty. 



For Kepler-8b, this approach was no possible since the 
r.m.s. scatter varies with time. To estimate the errors, we 
calculate the standard deviation using MAD in 50 point 
bins, excluding the transits and secondary eclipses. This 
binning window is shifted along by one data point and 
then repeated until we reach the end of the time series. 
The MAD-derived standard deviations are then plotted 
as a function of median time stamp from the binning 
window. Finally, we find that a quadratic fit through 
the data provides an excellent estimation of the trend. 
This formula is then used to generate all measurement 
uncertainties. 

3.3. Outliers 

Althou gh the data has bee n 'cleaned' by the Kepler 
pipeline ([Jenkins et al.ll20ioH ). we run the data through 
a highly robust outlie r detection algorithm developed by 
iBeaulieu et all (|2010l ) using the MAD method. Essen- 
tially, we use MAD to estimate the standard deviation 
and then estimate a sigma clipping level. This is found 
by setting the sigma clipping level to a value such that if 
all the errors were Gaussian, we would only expect to re- 
ject one valid data point by accident, which is a function 
of the array length. 

For each system we use MAD to estimate the outlier- 
robust standard deviation to be 100.7 ppm, 126.6 ppm, 
121.4 ppm and 96.4 ppm for Kepler-4b through 7b respec- 
tively. For Kepler-8b we find that the standard deviation 
of the photometry is variable as described in the previous 
subsection. We present the list of rejected outlier points 
in Table 2 of the appendix. 

3.4. Limb Darkening 

Accurate limb darkening coefficients were cal- 
culated for the Kepler bandpass for each planet. 
For the Kepler-bandpass, we used the high res- 
olution Kepler transmission function found at 
http://keplergo.arc.nasa.gov/CalibrationResponse.shtnill 
We adopted the SME-derived stellar properties reportea 
in the o rigina l discovery papers. We employed the 
iKurucd (|2006( ) atmosphere model database providing 
intensities at 17 emergent angles, which we interpolated 
linearly at the adopted Toff and logg values. The 
passband-convolvcd intensities at each of the e mergent 
angles were calculated following the procedure in iGlaretl 
(j2000( ). To compute the coefficients we considered the 
following expression: 



= l-^c,(l-/i)* 



(9) 



k=l 



where / is the intensity, fi is the cosine of the emer- 
gent angle, and Ck are the coefficients. The final coef- 
ficients resulted from a least squares singular value de- 
composition fit to 11 of the 17 available emergent angles. 
The reason to eliminate 6 of the angles is avoiding exces- 
sive weight on the stellar limb by using a uniform sam- 
pling (10 M values from 0.1 to 1 , plus /i — 0.05), as sug- 
gested bv iDiaz-Cordoves et al.l (|l995f ). This whole pro- 
cess is performed by a Fortran code written by I. Ribas 
(|Beaulieu et al.l [20101) . The coefficients are given in the 
final parameter tables for each planet. 
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3.5. Integration Time 

The Kepler long-cadence (LC) photometry is produced 
by the onboard stacking of 270 images of 6.02 seconds 
exposure. There is a 0.52 second read-time on-top of the 
exposure time which leads to a net duty-cycle of 91.4%. 
The duty cycle is sufficiently high that we can consider 
the LC photometry to be equivalent to one long exposure 
of 29.4244 minutes. As a result of such a long integration, 
sharp changes in the photometry, due to say a transit 
ingress, are smeared out into broader morphologies. In 
itself, this does not pose a problem since the nature of 
t his smear i ng is w ell understood and easily predicted. 

iKippind (|2010b( ) discusses in detail the consequences of 
finite integration times, with particular focus on Kepler's 
LC photometry. The effects of the smearing may be ac- 
counted for by using a time-integrated flux model for the 
lightcurve, as opposed to an instantaneous flux model. If 
we denote the flux predicted at any instant by F{t), the 

time- integrated flux by F{t) and the integration time by 
Z, the observed flux will be given by: 



Hi) 



ri+X/2 
Jt=i-X/2 



F{t)dt 



t+I/2 
t=i-I/2 



dt 



(10) 



IKippind ()2010bD proposes two methods of numerically 
integrating this expression; Simpson's composite rule and 
re-sampling. In this work, both of us compensate for the 
effect in different ways. Method A employs Simpson's 
composite rule and method B employ s re-sampling. I n 
both cases, we use the expressions of IKippind (I2010H) . 
given below in equations (5) and (6) to choose our nu- 
merical resolution such that the maximum photometric 
error induced by the finite numerical resolution is less 
than the observed photometric uncertainties, as seen in 
Table [II 



Comp.Simp. 
Resampling 



S I 
T 24m^ 



S I 
T 



(11) 

(12) 



, where 5 is the transit depth, r is ingress duration, and 
X is the integration time. 

4. KEPLER-4b 
4.1. Global Fits 



The d iscovery of Kepler-4b was made bv lBorucki et al.l 
(|2010bD . The planet is particularly interesting for joining 
the club of HAT-P-llb and GJ 436b as a Neptune-mass 
transiting exoplanet. Kepler-4b exhibits a sub-mmag 
transit around a 12**^ magnitude star, which leads to rel- 
atively large uncertainties on the system parameters but 
demonstrates the impressive performance of the Kepler 
photometry. 

However, the combination of the course sampling 
(e.g. the ingress/egress duration is 3-4 times shorter than 
the cadence of the observations) , the very low RV ampli- 
tude and a sub-mmag transit make Kepler-4b the most 
challenging object to determine reliable system parame- 
ters for in this paper. The BIC test for the approriate 
RV model prefers the simplest description possible, re- 



flecting the low signal-to-noise levels encountered for the 
radial velocities. 

The largest di f ference between our own fits and that 
of iBorucki et al.l (l2010bf) is the retrieved a/i?* and im- 
pact parameter, h. IBorucki "eTall (|2010tf) find an equa- 
torial transit solution whereas our method A circular fit, 
and both modes of method B, place Kepler-4b at an im- 
pact parameter of 0.5-0.6. Curiously, the eccentric fit of 
method A prefers an even larger impact parameter than 
this, as a result of letting the limb darkening be fitted. 

Due to t he well-known neg ative correlation between b 
and a/Ri, (jCarter et al.ll2008l ). these larger b values lead 
to lower a/Ri, values and consequently a significantly 
lower lightcurve derived stellar density, p*. Indeed, the 
lower stellar density results in one of the largest stellar 
radii out of all the known transiting systems. The fitted 
models on the phase-folded lightcurves are shown on Fig- 
ure [2] using method B. Correlated noise was checked for 
in the residuals using the Durbin- Watson statistic, which 
finds d — 2.082, well inside the 1% critical boundary of 
2.135. The orbital fits to the RV points are shown in 
Fig. [3] (method B), depicting both the circular and ec- 
centric fits. The final table of all results are shown in 
Tabled 

4.1.1. Eccentricity 

From Table [T4l of the Appendix, the circular fit with no 
other free parameters is the preferred model to describe 
the radial velocities of Kepler-4b. We may perform other 
tests aside from the BIC model selection though. The 
global circular fit, using method A, yields a — 28.54 
and the eccentric fit obtains — 20.56. These values 
correspond to the lowest solution of the simplex global 
fit, but for the RV points only, which dominate the eccen- 
tricity determination. Assuming the period and epoch 
are essentially completely driven by the photometry, the 
number of free parameters are 2 and 4 for the circular 
and eccentric fits respectively. Therefore, based upon 
an F-test, the inclusion of two new degrees of freedom 
to describe the eccentricity is justified at the 1.7-cr level 
(91. 5% confidence). Anoth er test we can implement is 
the iLucv k, SweenevI ()197lD test, where the significance 
of the eccentric fit is given by: 



P(e > 0) = 1 - exp 



2 n 



e 
2&i 



(13) 



Where e is the modal value of e and cr^ is the error 
(in the negative direction). We find that this gives a sig- 
nificance of 88.2% or 1.6-(T, in close agreement with the 
F-test result. We also computed the posterior distribu- 
tion of the distance of the pericentre passage, in units 
of the stellar radius, and found (a/i?*)(l — e) = 2.6lo'g 
for method A. Therefore, if the eccentric fit is the true 
solution, Kepler-4b would make an extremely close peri- 
centre passage. 

We note that IBorucki et al.l (|2010b[ ) found a best-fit 
eccentricity of e = 0.22 (no quoted uncertainty) but the 
authors concluded the eccentric fit was not statistically 
significant, but provided no quantification. Based upon 
the current observations, it is not possible to make a 
reliable conclusion as to whether Kepler-4b is genuinely 
eccentric or not. At best, it can be considered a ~ 2- 
(7 marginal detection of eccentricity. The only assured 
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TABLE 1 

Integration resolutions required to match observed photometric performance, based upon expressions of \Kivvin3l f201 OH) . 



Planet 


Required m for 


m used for 


Required N 


N used 




Comp. Simp. Rule 


Comp. Simp. Rule 


for Resampling 


for Resampling 


Kepler-4b 


0.9434 


2 


1.6339 


4 


Kepler-5b 


1.7108 


2 


2.9632 


4 


Kepler-6b 


2.1848 


2 


3.7839 


4 


Kepler- 7b 


1.7199 


2 


2.9790 


4 


Kepler-8b 


1.0810 


2 


1.8723 


4 



thing we can say is that e < 0.43 to 95% confidence. 
Despite the ambiguity of the eccentricity, it is interesting 
to consider the consequences of this object possessing 
e > 0. 

For the eccentric fit, the extreme proximity of this 
planet to the star raises the issue of tidal circulariza- 
tion. Let us continue under the assumption that no 
third body is responsible for pumping the eccentricity 
of Kepler-4b. For a planet initially with e '--^ 1, the ec- 
centricity decreases to 1/ exp(l)" after n circularization 
timescales. This means the maximum number of circu- 
larization timescales which have transpired is given by 
n < — log(e). Since n = Tage/rcirc (age of the system 
divided by circularization timescale) then we may write 
Tcirc ^ Tagc/n. Usiug thc mcthod of I Adams fc LaughlinI 
(|2006i) . the tidal circularization time may be written as: 



2\13/2 



(14) 



Given that Tcirc ^ Tg^gc/n, we may compute the mini- 
mum possible value of Qp through re-arrangement; 



Qp> 



63 2tt 



'Ip\{a/R^)J 



'age 



1 



(15) 



log(e) (1 -e2)i3/2 

The final term, (1 — e^)^'^/^ may be neglected since 
we are interested in the minimum limit of Qp and this 
term only acts to further increase the Qp for non-zero 
e. The advantage of doing this is that e is a function 
of time and thus we can eliminate a term which would 
otherwise have to be integrated over time. Taking the 
posterior distribution of the Q p equation for method A 
yields Qp > (Lll^g) x 10^. We find that 99.9% of the 
MCMC trials correspond to Qp > 10^ for the eccen- 
tric fit. Thus, if the eccentric fit was genuine, Kepler- 
4b would exhibit very poor tidal dissipation. This is 
somewhat remin iscent of G J 436b ()Deming et al.l [20071) 
and HAT-P-llb ()Bakos et al.ll2009[ ) for which a large Qp 
value is also predicted and raises the interesting question 
as to whether hot-Neptune generally possess larger Qp 
values than their Jovian counterparts. We also note that 
WASP- 18b has a lso bee n proposed to possess a greater- 
than-Jovian-Qp (jHellier et al- 2009 ). 

This discussion highlights the importance of obtaining 
more statistics regarding the eccentricity of transiting 
hot-Neptunes, which may reveal some fundamental in- 
sights into the origin of these objects. 

4.1.2. Secondary eclipse 

We here consider the circular model only since the ec- 
centric fit is not unambiguously accepted. Global fits 



do not find any significant detection of the secondary 
eclipse. Secondary depths > 104 ppm are excluded to 
3-cr confidence, which corresponds to geometric albedos 
greater than unity and thus this places no constraints 
on the detectability of reflected light from the planet. 
The thermal emission is limited by this constraint such 
that Tp^brightness < 3988 K, which places no constraints 
on redistribution. 



1.000 



v: 



0.999 



g. 0.998 



0.997 



.'.•^i'mtII' 




l|'!|HlitH|i|i!*jt|il|iH|ii*l}iM|^^^ 



-0.04 -0,02 0,00 0,02 0,04 

Orbital phsse 

Fig. 2. — Phase-folded primary transit lightcurve of Kepler-4 
for method B. The upper curve shows the eircular fit, the bottom 
eurve the eecentric fit. Solid (blue) lines indieate the best fit resam- 
pled model (with bin- number 4). The dashed (red) lines show the 
eorresponding unbinned model, which one would get if the transit 
was observed with infinitely fine time resolution. Residuals from 
the best fit resampled model for the circular and eccentric solutions 
are shown below in respective order. 



4.1.3. Properties of the parent star Kepler-4 

The Yonsei-Yale model isochrones from IYi et al.l (|2001[ ) 
for metallicity [Fe/H]=-|-0.17 are plotted in Figure 21 
with the final choice of effective temperature T^e* and 
a/i?i marked, and encircled by the la and 2a confidence 
ellipsoids, both for the circular and the eccentric orbital 
solution from method B. Here the MCMC distribution of 
a/Ri, comes from the global modeling of the data. Nat- 
urally, errors of the stellar parameters for the eccentric 
solution are larger, due to the larger error in a/i?^, which, 
in turn, is due to the uncertainty in the Lagrangian or- 
bital parameters k and h (as opposed to fixing these to 
zero in the circular solution). 

The stellar evolution modeling provides color indices 
that may be compared against the measured values as 
a sanity check. The best available measurements are 
the ne ar-infrared magnitud es from the 2MASS Cata- 
logue ijSkrutskie et al.| [2006[ ). For Kepler-4, these are: 
^2MAss = 11-122 ± 0.023, iJ2MASS = 10.836 ± 0.031 and 
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20 




Phase with respect to T^, 

Fig. 3. — Top: RV measurements from Keck for Kepler-4, along 
with an orbital fit, shown as a function of orbital phase, using 
our best-fit period. Solid (blue) line shows the circular orbital fit 
with binned RV model (3 bins, separated by 600 seconds). The 
dashed (red) line shows the eccentric orbital fit (with the same 
bin parameters). Zero phase is defined by the transit midpoint. 
The center-of-mass velocity has been subtracted. Note that the er- 
ror bars include the stellar jitter (taken for the circular solution), 
added in quadrature to the formal errors given by the spectrum 
reduction pipeline. Fits from method B. Middle: phased residu- 
als after subtracting the orbital fit for the circular solution. The 
r.m.s. variation of the residuals is 3.8ms~^, and the stellar jitter 
that was added in quadrature to the formal RV errors is 2.0 ms~^. 
Bottom: phased residuals after subtracting the orbital fit for the 
eccentric solution. Here the r.m.s. variation of the residuals is 
3.5ms^^, and the stellar jitter is 1.6ms~^. 




6000 5500 5000 

Effective temperature [K] 

Fig. 4. — Stellar isochrones from IYi et al.l 11200 ll *) for metallicity 
[Fe/H]=-(-0.17 and ages 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0, 7.5, 8.0, 
8.5 and 9.0 Gyrs. The final choice of T^ffi, and a/Ri, for the circular 
solution are marked by a filled circle, and is encircled by the Icr and 
2(T confidence ellipsoids (solid, blue lines). Corresponding values 
and confidence ellipsoids for the eccentric solution are shown with 
a triangle and dashed (red) lines. Fits from method B. 

.f^2MASS = 10.805 ± 0.023; which we have converted to 
the photometric system of the model s (ESQ system) us- 
ing the transformations bv lCarpenteii ()2001[ ). The result- 
ing measured color index is J — = 0.340 ± 0.009. This 
is within 2-a of the predicted value from the isochrones 
oi J ~K = 0.38 ± 0.02 (seeHD. The distance to the ob- 
ject may be computed from the absolute K magnitude 
from the models (Mk = 1.95 ± 0.35) and the 2MASS 
Ks magnitude, which has the advantage of being less af- 
fected by extinction than optical magnitudes. The result 



is 600lg5 pc, where the uncertainty excludes possible 
systematics in the model isochrones that are difficult to 
quantify. 

4.2. Transit timing analysis 
4.2.1. Analysis of variance 

We find TTV residuals of 207.8 seconds and TDV 
residuals of 217.0 seconds using method A. Repeating the 
TTV analysis for method B finds best-fit times of con- 
sistently less than 0.25-ct deviance (we note that similar 
results were found for all planets). The TTV indicates 
timings consistent with a linear ephemeris, producing a 
of 5.7 for 11 degrees of freedom. This is somewhat 
low with a probability of 10.9% of occurring by coin- 
cidence. This may be an indication that the MCMC 
methods adopted here yield artificially large error bars, 
perhaps due to our choice of free-fitting every lightcurve 
rather than fixing certain parameters. However, we pre- 
fer to remain on the side prudence by overestimating such 
errors rather than underestimating them. The TDV too 
is consistent with a non-variable system exhibiting a 
of 6.4 for 12 degrees of freedom (10.6% probability of 
random occurrence). 

Under the assumption that the timing uncertainties 
are accurate or overestimated, it is possible for us to 
determine what signal amplitudes are excluded for wave- 
forms of a period less than the temporal baseline. For 
the TTV, we note that for 11 degrees of freedom, excess 
scatter producing a of 28.5 would be detected to 3-cr 
confidence. This therefore excludes a short-period signal 
of r.m.s. amplitude > 341.6 seconds. Similarly, for the 
TDV, we exclude scatter producing a x^ of 30.1, or a 
short-period r.m.s. amplitude of > 339.3 seconds, to 3- 
(7 confidence. This constitutes a 4.9% variation in the 
transit duration. 

These limits place constraints on the presence of per- 
turbing planets, moons and Trojans. For the case a 2:1 
mean- motion resonance perturbing planet, the libration 
period would be ~ 328 cycles (2.9 years) and thus we do 
not possess sufficient baseline to look for such perturbers 
yet. However, this libration period is less than the mis- 
sion lifetime of Kepler (3.5 years) and so assuming the 
same timing precision for all measurements, the TTV 
will be sensitive to a 0.14 perturber. 

The current TTV data rules out an exomoon at the 
maximum orbital radius (for retrograde revolution) of 
mass > 11.OM0. The TDV data rules an exomoon of 
> 13.3 M0 at a minimum o rbital distance. 

Using the expressions of IFord fc HohIIa3 (pOOl and 
assuming Mp ^ Mxrojan, the expected libration period 
of Kepler-4b induced by Trojans at L4/L5 is 50.5 cycles 
and therefore we do not yet possess sufficient temporal 
baseline to look for the TTVs. However, inspection of 
the folded photometry at ±P/6 from the transit centre 
excludes Trojans of effective radius > 1.22i?0 to 3-a con- 
fidence. 

4.2.2. Periodograms 

As discussed in section §2.3.2, we may use an F-test pe- 
riodogram to search for periodic signals, as this negates 
any problems introduced by underestimating or overes- 
timating timing uncertainties. Whilst a search for ex- 
cess variance is sensitive to any period below the tem- 
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TABLE 2 

Global fits for Kepler-4b. Quoted values are medians of MCMC trials with errors given by 1-a quantiles. Two independent methods are 
used to fit the data (A&B) with both circular and eccentric modes (c&e). * = fixed parameter; f = parameter was floated but not fitted. 



Parameter 


Discovery 


Method A.c 


Method A.e 


Model B.c 


Model B.e 


Fitted params. 


P/days 
E (BJD-2,454,000) 


3.21346 ± 0.00022 
956.6127 ± 0.0015 


O.zio^0_o 00058 

956.61251HS11 


n Qi Q^Q+U.UUU3U 
O.ziot0_o 00030 

956.61322l°:0027 


3.21344 + 0.00024 
956.6124 ± 0.0017 


3.21340 + 0.00027 

956.6132 + 0.0020 


Ti,4/s 




142591^°! 


142221^35 


14282 + 190 


14187 + 207 


T/s 




135781593 


13550l33fo 


13733 + 182 


13706 + 204 


(Tl,2 ~ T3,4)/S 




c:io+614 
J^>-|_149 


80911°° 


501 + 207 


475 + 104 


62 

(C/iJ*)/days-i 




0691+° °"^ 
U.UDyi_o 0070 

0.30l°:i 
12.881°:!? 


0740+° °°94 
U.Ur4:U_o 0063 
54+0-19 

ici.OU_o 5g 


0.0660 + 0.0046 

o.29oiE;:i«f, 

12.58 + 0.17 


0.0656 + 0.0031 

0.25910:^32 

12.61 + 0.20 


(i^P,d/-F*)/ppm 
esinw 
ecoso) 


0* 
0* 

■-'•''-1.9 


10+29 
1U_29 

0* 
0* 

9.2l3;3 


0+67 
^-355 
24+°- 12 

033+° °5* 

U.UOO_o 061 


-6 + 13 
0* 
0* 

9.3 + 1.4 


5 + 29 
0.06210:11 
0.083lo:j3g 
9.1 + 1.5 




-1.27± 1.1 


-0.9t2;2 


1-^-1.1 


-0.9 + 0.9 


-0.8 + 1.0 


'y/ms~lday~l 


0* 


0* 


0* 


0* 


0* 


ttroj /days 


0* 


0* 


0* 


0* 


0* 


B 


1.02 ±0.02 


1.02 ± 0.02 t 


1.02 + 0.02 t 


1.02 + 0.02 t 


1.02 + 0.02 t 


SME derived params. 


Teff/K 


5857 ± 120 


5857 ± 120 


5857 + 120 


5857 + 120 


5857 + 120 


log(g/cgs) 


4.25 ±0.10 


4.25 ±0.10 


4.25 + 0.10 


4.25 + 0.10 


4.25 + 0.10 


[Fe/H] (dex) 


+0.17 ±0.06 


+0.17 + 0.06 


+0.17 + 0.06 


+0.17 + 0.06 


+0.17 + 0.06 


II sin j (kms~l) 


2.2 ±1.0 


2.2 + 1.0 


2.2 + 1.0 


2.2 + 1.0 


2.2 + 1.0 


Wl 

U2 




61 +"-59 
U.01_o 39 

-0 91+0-52 
"•^-'■-0.68 


n t;4+0.63 

U.04_o 35 

-0 23+° ''^ 


0.4086* 
0.2633* 


0.4086* 
0.2633* 


Model indep. params. 


Hp/ R-k 
a/Ri, 
b 

il° 


u.uz4:ru_Q 00030 

fj 47+0.26 
"•^'-0.28 

0.022l°;234 


0263+" "°22 
U.UZDd_Q 0014 
K 4Ci+0.87 

o.55l°;f5 

84.3t^;° 


0272+° °°l ' 
u.uz(z_o 0012 

„+1.34 

'^■'-"-0.63 

0.73l°i2 
74.ll?-5 


0.0257 + 0.0009 
5 41+0.56 

O-^l_0.82 
5Qq+0.155 
u.dO»_o 205 

84.31^:^ 


0.0256 + 0.0006 

c -,7+1.50 
^■-^ ' -0.94 

0.50910:113 

84.2111 


e 


0* 


0* 




0* 


0.208 + 0.104 








84.5111:? 




76 + 145 


RV jitter (ms-i) 
p*/gcm-3 

log(3p/cgs) 
5'i,4/s 
tsec/(B.ID-2, 154,000) 


0.525 

o 1(5+0.06 


1.7 

'^■■"J-0.18 

2-90lo.38 
14259t™ 

961.4325l°:°°«« 


0.08l+0:J34 
2.58l°:^« 

r,+14946 
'^0 

961.28l?g? 


2.0 
0.29 + 0.10 

2 96+0-13 

14282 + 190 
980.713 + 0.001 


1.6 
0-251°:?? 
2.92 + 0.19 
11975 + 544 
980.89 + 0.17 


Derived stellar params. 


R*/Rq 
log(5/cgs) 

My /mag 
Age/Gyr 

Distancc/pc 


1 r,Qn+0.U53 

l.zzo_o 091 

l-4='-0.()84 

4.17 ±0.04 

O ofi+0.66 
^•^"-0.48 

4.00 ± 0.28 
4.5 ±1.5 
550 ± 80 


l-^°-0.13 

1 09+0. 81 
l-°^-0.28 
4 02+"-12 

3.5+?-? 

4t;+0.40 

49+2.1 
^■■'-1.2 

566l2f 


1 56+0-20 
l.dO_o.24 

2 07+0.79 

-0.93 

3.6810-25 

q q+5.7 
y.o_4.9 
n QQ+0.82 
z.oa_Q 52 

2 74+1-31 
z.'^_0.86 
Q99+252 
»^-'_289 


1.275 + 0.109 

1 S9Q+0-426 

4.01 + 0.11 

3 56+2-03 
"-'■"-'"-0.77 

3.43 + 0.37 

4.6li:J 
eool^f 


1 909+0. 160 
l.zoz_Q -1^23 
1 Q99+0-572 
i.yzz_Q 430 

3.98 + 0.18 

3.9il?:f2 

3.33 + 0.58 

4 s+i-8 
^■"^-1.0 

630ll|« 


Derived planetary params. 


Mp/M} 
Rp/Rj 


0.077 ±0.012 
0.357 ±0.019 

qi+0.36 


U.U»l_o 030 
U.'iDU_o 084 


0.09610:02^ 
0.7910:2? 

0.24l°:f3 


0.079 + 0.013 
n 4^7+0.134 

U.^Or_o 056 

1.00 + 0.49 


0.077 + 0.015 
0.480l0:l« 
0-87lJ:i| 


a/AU 


0.0456 ± 0.0009 


0.0463tO:SE'23 


n 0494+0-0020 
U.U4y4_o 0026 


0.0462 + 0.0013 


0.0463 + 0.0017 


Teq/X 


1650 ± 200 


1777l30« 


991 c: + 233 
zzid_339 


1783l^f 


1832 + 212 
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poral basehne, periodograms are hmited to a range of 

'^Pp ^ -^signal ^ 2_f\)ascliiic- 

There are no peaks in the TTV periodogram below 
5% FAP. The TDV does seem to exhibit a signifi- 
cant peak, which can also be seen by eye in the TDV 
plot itself, at a period of (15.5 ± 3.1) cycles, ampli- 
tude (204.9 ± 61.2) seconds and FAP 2.1% (2.3-cr). A 
linear trend in the durations could be confused with 
a long period signal. Fitting a simple a + bt model 
yields a = (170.9 ± 95.2) seconds and b ^ (-31.6 ± 
13.8) seconds/cycle with FAP 2.6% (2.2-cr). Therefore, 
the two models yield very similar significances and it 
is difficult to distinguish between these two scenarios. 
The linear drift model seems improbable due to the very 
rapid change in duration the trend implies. The gra- 
dient equates to a change of around one hour per year 
in the half-duration, i.e. two hours per year in T. The 
period of the signal is too large to be caused by an ex- 
omoon and TDV is generally insensitive to perturbing 
planets, meaning that any planet which could produce 
such a large effect would dominate the radial velocity. 
Other TDV effects, for example apsidal precession, are 
not expected to occur on timescales of 10-20 days which 
puts the reality of the signal in question. 

The phasing of the measurements, as defined in §2.3.4, 
yields a clear periodic waveform of period 4 cycles. 
Therefore one possible explanation for the TDV signal 
is a mixing of this phasing frequency with the Nyquist 
frequency to produce a signal of period 16 cycles, which is 
consistent with the value determined earlier. The analy- 
sis of further measurements is evidently required to assess 
the reality of the putative signal. 

4.3. Depth and OOT Variations 

The transit depth is consistent with the globally fitted 
value yielding a of 4.3 for 12 degrees of freedom. Sim- 
ilarly, the OOT levels are flat yielding a of 6.1 for 12 
degrees of freedom. 

5. KEPLER-5b 

5.1. Global fits 

5.1.1. Comparison to discovery paper 

Kepler-5b was discovered by iKoch et al.l (|2010[ ) . Our 
global fits reveal Kepler-5b to be a ^ 1.35i?j planet on 
an orbit consistent with that of a circular orbit, tran- 
siting the host star at a near-equatorial impact parame- 
ter. The fitted models for method B of the phase-folded 
lightcurves are shown on Figure [Bl Correlated noise was 
checked for in the residuals using the Durbin- Watson 
statistic, which finds d = 1.991, which is well within 
the 1% critical level of 1.867. The orbital fits to the RV 
points are shown in Fig. |9] for method B, depicting both 
the circular and eccentric solutions. We obtain transit 
parameters broadly consistent with that of the discovery 
paper, except for b where method A favors an equato- 
rial transit and method B favo rs a near-equatorial tran- 
sit, whereas iKoch et al.l (|2010( ) found 6 ~ 0.4. The final 
parameters are summarized in Table |4l Lightcurve fits 
reveal that the theoretical limb darkening values differ 
from the fitted values by a noticeable amount and the 
residuals of method B do reveal some structure as a con- 
sequence. Future SC mode observations will pin down 
the limb darkening to a much higher precision. 



5.1.2. Eccentricity 

The global circular fit of method A yields a = 20.77 
for the RVs and the eccentric fit obtains = 19.77. 
Based on an F-test, the inclusion of two new degrees of 
freedom to describe the eccentricity is justified below the 
l-fj level and is therefore not significant. We therefore 
conclude that the system is consistent with a circular 
orbit and constrain e < 0.086 to 95% confidence. 

5.1.3. Secondary eclipse 

We here only consider the circular fits since the ec- 
centric solution is not statistically significant. Method 
B provides the more accurate constraints due the use of 
multiple baseline scalings. Proceeding with method B, 
the secondary eclipse is weakly detected in our analysis 
to 93.4% confidence, or 1.8-cr, with a depth of 26ti[7ppm. 

The MCMC distribution of the Fp^d/F* values is shown 
in Figure [51 The discovery paper reported a weak 2.0-tT 
detection of the eclipse but did not provide an amplitude 
for the signal. An eclipse of depth > 76 ppm is excluded 
at the 3-cr level, which excludes a geometric albedo > 
0.47 to the same level. 

If due to thermal emission, we find that by integrat- 
ing a Planck function with the Kepler response function, 
the planetary brightness temperature would have to be 
2480^284 which is inconsistent with that expected for 
the equilibrium temperature of (1770 ± 19) K. Assum- 
ing negligible albedo and solving for the redistribution 
factor implies a factor of 3.7, which surpasses the maxi- 
mal physically permitted value of 8/3 (corresponding to 
2300K). Internal heating also seems to be unlikely as the 
heat source would have to be able to generate a sur- 
face temperature of 2300 K alone, or have an intrinsic 
luminosity of 4.6 x 10~'*Lq which is approximately that 
of a M7V star. Tidal heating from eccentricity seems 
improbable given that the orbit is consistent with a cir- 
cular orbit. Ho wever, using our best-f it eccentricity and 
the equations of lPeale fc CassenI ()1978[ ). we estimate that 
tidal heating could induce 8.3 x 1O~^L0 for Qp = 10^ 
and = 0.5. In order to get enough tidal heating, 
we require Q' = Qpjliip = 35, which is much less than 
the typical values expected. Furthermore, such a low 
Qp leads to very rapid circularization which somewhat 
nullifies the argument. 

If the eclipse was due to reflection, we estimate a re- 
quired geometric albedo of Ag = 0.15 ± 0.10. This latter 
option does offer a plausible physical origin for the sec- 
ondary eclipse and therefore should be considered as a 
possible explanation for the observation. Although there 
has been a general impression that hot-Jupiters must 
have low albedos after the remarkable MOST constraints 
oi Ag < 0.17 for HD 209458b by R owe et all (|2008f ). a 
recent study by Cow an fc Agol (20ld ) finds evidence for 
hot-Jupiters having much higher albedos (> 0.4) in a sta- 
tistical sample of 20 exoplanets. Therefore, a geometric 
albedo of 0.15 is not exceptional in such a sample and 
offers much more reasonable explanation that thermal 
emission. However, at this stage, the significance 
is too low to meet the requirements for formal detection. 

5.1.4. Properties of the parent star Kepler-5 

The Yonsei-Yale model isochrones from iYi et al.l (|2001l) 
for metallicity [Fe/H]=-l-0.04 are plotted in Figure [TOl 
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Fig. 5. — Upper Left: TTV (squares) and TDV (diamonds) for Kepler-J^b (see ^2.3.1 for details). Upper Right: TTV periodogram 
(solid) and TDV periodogram (dashed) for Kepler-4-b, calculated using a sequential F-test (see ^2.3.2 for description) . Lower Left: 
Transit depths from individual fits of Kepler-4b (see ^^2.3.1 for details). Lower Right: "Phasing" of Kepler long cadence photometry for 
Kepler-4b (see ^2.3.4 for description). 



TABLE 3 

Mid-transit times, transit durations, transit depths and out-of-transit (normalized) fluxes for Kepler-4b. 



Epoch tc/{BJD-2,454,000) T/s (R^jRlfJX F\^t 
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Fig. 6. — Top: Phase-folded primary transit lightcurve of 
Kepler-5. The upper curve shows the circular fit, the bottom curve 
the eccentric fit. Sohd (blue) lines indicate the best fit resampled 
model (with bin- number 4). The dashed (red) lines show the corre- 
sponding unbinned model, which one would get if the transit was 
observed with infinitely fine time resolution. Residuals from the 
best fit resampled model for the circular and eccentric solutions 
are shown below in respective order. 
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Fig. 7. — Phase-folded secondary transit lightcurve of Kepler-5 
(top), and residuals from the best fit (bottom). Only the fit for 
the circular orbital solution of method B is shown. 




2e-05 4e-05 6e-05 

Planet dayside flux / stellar flux 

Fig. 8. — Distribution of Fp^^/Fi, for Kepler-5from the global 
modeling. Fp^d/Ft, is primarily constrained by the depth of the 
secondary eclipse. Only results for the circular orbital solution of 
method IS are shown. 

with the final choice of effective temperature Toff* and 
a/Ri, marked, and encircled by the la and 2a confidence 
ellipsoids, both for the circular and the eccentric orbital 
solution. Here the MCMC distribution of a/i?* comes 
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Phase with respect to T^, 

Fig. 9. — Top: RV measurements from Keck for Kepler-5, along 
with an orbital fit, shown as a function of orbital phase, using 
our best-fit period. Solid (blue) line shows the circular orbital fit 
with binned RV model (3 bins, separated by 600 seconds). The 
dashed (red) line shows the eccentric orbital fit (with the same 
bin parameters). Zero phase is defined by the transit midpoint. 
The center-of-mass velocity has been subtracted. Note that the er- 
ror bars include the stellar jitter (taken for the circular solution), 
added in quadrature to the formal errors given by the spectrum 
reduction pipeline. Fits from method B. Middle: phased residu- 
als after subtracting the orbital fit for the circular solution. The 
r.m.s. variation of the residuals is 16.8 ms~^, and the stellar jitter 
that was added in quadrature to the formal RV errors is 15.5 m s~^. 
Bottom: phased residuals after subtracting the orbital fit for the 
eccentric solution. Here the r.m.s. variation of the residuals is 
16.3 ms^^, and the stellar jitter is 15.0 ms~^. 

from the global modeling of the data. 
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Fig. 10. — Stellar isochrones from lYi et al.l l|20011 ) for motal- 
licity [Fo/H]=-f0.04 and ages 1.0, 1.4, 1.8, 2.2, 2.6, 3.0, 3.4 and 
3.8Gyrs. The final choice of T^ff^. and a/i?* for the circular solu- 
tion are marked by a filled circle, and is encircled by the Icr and 
2(T confidence ellipsoids (solid, blue lines). Corresponding values 
and confidence ellipsoids for the eccentric solution are shown with 
a triangle and dashed (red) lines. Fits from method B. 



5.2. Transit timing analysis 
5.2.1. Analysis of variance 

We find TTV residuals of 30.1 seconds and TDV resid- 
uals of 35.6 seconds. The TTV indicates timings consis- 
tent with a linear ephemeris, producing a of 8.7 for 10 
degrees of freedom. The TDV is also consistent with a 
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TABLE 4 

Global fits for Kepler-5b. Quoted values are medians of MCMC trials with errors given by 1-a quantiles. Two independent methods are 
used to fit the data (A&B) with both circular and eccentric modes (c&e).* = fixed parameter; f = parameter was floated but not fitted. 



Parameter 



Discovery 



Method Ax 



Method A.e 



Model B.c 



Model B.e 



Fitted params. 
P/days 
E (BJD-2,454,000) 

T/s 

(Tl,2 ~ T3,4)/S 

(C/i?*)/days-i 
(Fp^d/F*)/ppm 

esinoj 

ecoso) 

7/ms~-' 

'y/ms^-'-day"-'- 
ttroj /days 
B 



955.90126l°;rof9 

1+184 
-165 

15209liS 
1-257 
-84 
7-1-0.026 
'-0.017 
^4-0.143 



955.90124l°:°~ii 

15232tJig 
127511^3 

-1-0.018 
-0.012 
-1-0.106 
0.038 



3.548460 ± 0.000032 
955.90122 ± 0.00021 



16580Z 
15209; 
1309j 
0.637"'' 



0.066] 



().or,9 

, + 0.14 



ll-42lo.l3 

-til 



0.632 
0.042 

11 403+«-09» 



29Z 



0* 
0* 

231.3 ±2.8 
-46.7 ± 4.1 

0* 

0* 

1.020 ±0.002 



0* 

0* 

220.6ti4;2 



-43.1"' 



hio.o 

-10.0 
28+°-2® 



+24 
24 
+0.041 
-0.038 
+0.017 
0.016 
221.ll|i 
-42.1+«-» 



29 
0.012 
-0.008 



0.28 



7.0 
+0.17 
0.17 



0* 

1.020 ± 0.002 t 



0* 

1.020 ± 0.002 t 



3.548459 ± 0.000035 
955.90124 ± 0.00023 
16459 ± 35 

15233 ± 23 
1227 ± 26 
0.6225 ± 0.0032 

11.34 ±0.02 
26± 17 
0* 
0* 

228.1 ±8.1 
-44.8 ± 6.5 

0* 

0* 

1.020 ±0.002 t 



3.548459 ± 0.000035 

955.90124 ± 0.00023 

16468 ± 35 

15230 ± 25 

1227 ± 35 

0.6225 ± 0.0032 

022+° °3° 
"■"^^-0.013 

11.35 ±0.02 

25 ± 25 
0.008 ± 0.044 
-0.009 ± 0.023 
228.5 ±8.1 
-43.9 ± 6.8 
0* 
0* 

1.020 ± 0.002 t 



SME derived params. 



TeB/K 


6297 ± 60 


6297 ± 60 


6297 ± 60 


6297 ± 60 


6297 ± 60 


log(g/cgs) 


3.96 ±0.10 


3.96 ±0.10 


3.96 ±0.10 


3.96 ±0.10 


3.96 ±0.10 


[Fe/H] (dex) 


+0.04 ± 0.06 


+0.04 ± 0.06 


+0.04 ± 0.06 


+0.04 ± 0.06 


+0.04 ± 0.06 


vsini (kms~^) 


4.8 ±1.0 


4.8 ± 1.0 


4.8 ±1.0 


4.8 ±1.0 


4.8 ± 1.0 


Ml 
U2 




"•^'^-0.12 
O-+0.25 
"••^'-0.27 


0.244l0;»«t 
0.39iri| 


0.3462* 
0.2972* 


0.3462* 
0.2972* 


Model indep. params. 


Ftpl Fi* 
a/R* 
b 

ir 


U.U0190_Q 00047 

6.06 ±0.14 
n QQCj+0.051 

86.3 ±0.5 


o.o798i:j;ijr,i? 

O.OOl'^;^^ 


U.U/yOZ_0 00074 
6 1ci+0-30 
D-lO_0.35 

^•"Jl-0.01 

SS.lll'l 


0.0789 ± 0.0002 

(, 1 40+0.067 
U.i4J_o 057 

88.7 ±0.6 


0.0789 ± 0.0002 
6.26 ± 0.29 
u.iia_o 060 


e 


0* 


0' 


U.UCS4_0 018 


0* 


0.038 ±0.028 


u;/° 






132lif 




150 ± 83 


RV jitter (ms"!) 


0.336 


8 2 
U.JDU_o 


7.7 

U.OOU„o 057 


15.5 

0.38t°:°^ 


15.0 

n q7+0.06 
"•■^'-0.04 


log(9p/cgs) 


3.41 ± 0.03 


q .qq+0.046 

0.'±00_Q 073 


q ^qq+0.048 

o.^oo_o og2 


3.48 ±0.02 


3.47 ± 0.04 


Si,4/s 




165801^8* 


20006ll|i9 


16459 ± 35 


16718 ± 1417 


tsec/(BJD-2,454,000) 




qfji 22^72+° °°°"^ 


959.43tJ;JI 


978.9660 ± 0.0002 


975.397 ± 0.052 



Derived stellar params. 



M*/Mq 


1 q-.+O.IMO 

l.O('±_o 059 


l.daj_o 03B 


1 q7r,+U.U48 

l-^'^-0.041 


1.352 ±0.026 


1.358 ±0.039 


R*/Rq 


1 7Qq+0.043 
1. /yo_o 052 


'^^-0.059 


^■'0»-0.094 


1 709+0.029 
1. HJ»_o 017 


1.730 ±0.091 


log{9/cgs) 


4.07 ± 0.02 


4-U87_o.o57 


4.USU_o.o48 


4.10 ±0.01 


4.09 ±0.04 


L*/Lq 


4 67+0-63 
-0.59 


4 34+0.77 
^•'-'^-0.37 


4 41+0.69 


4.12 ±0.22 


4.21 ±0.48 


Mv 1 mag 




o 1«+0.10 
•^•^^"-0.18 


, 14+0.14 
•^•-'-^-0.16 


3.22 ±0.07 


3.20 ±0.13 


Age/Gyr 


3.0 ±0.6 


Q oo~l~0.23 
z.oo_Q 25 


r, 09 + 0.24 
^•»2_o 24 


2.9 ±0.2 


2.9 ±0.2 


Distancc/pc 




iioil^t 


1109lg7 


906 ± 17 


916 ± 49 



Derived planetary params. 



^■-'--'-^-0.059 
1 401+0.041 
l.^Ji_0 052 

0.894 ± 0.079 
0.05064 ± 0.00070 
1868 ± 284 



r+U.i4 

3-0.14 
1 qr:9+0.149 
l.ooz_o 052 

1 01+°-^^ 
l->Jl_o.25 

0.05056tO;»?»^i 



9 Qr4+0.U99 
z.uoi_o.098 
^+0.112 
-0.076 
1+0.18 
-0.20 



Mp/Mj 

pp/gcm-^ 

a/AU 



2.05"* 



1790"* 



1.366: 
1.00"f 



u.uouoy_o. 00051 



1796 



+ 57 
45 



2.094 ±0.079 

1 qi 9+0.026 
1.O1Z_0 014 

1.15 ±0.06 
0.0503 ± 0.0003 
1770 ± 19 



2.101 ±0.087 
1.330 ±0.071 

1 11+0.19 
1-11-0.14 

0.0504 ± 0.0005 
1779 ± 43 
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non- variable system exhibiting a of 10-5 for 11 degrees 
of freedom. 

For the TTV, we note that for 10 degrees of freedom, 
excess scatter producing a of 26.9 would be detected 
to 3-CT confidence. This therefore excludes a short-period 
signal of r.m.s. amplitude > 39.8 seconds. Similarly, for 
the TDV, we exclude scatter producing a of 28.5, or 
a short-period r.m.s. amplitude of > 45.9 seconds, to 3- 
a confidence. This constitutes a 0.60% variation in the 
transit duration. 

These limits place constraints on the presence of per- 
turbing planets, moons and Trojans. For the case a 2:1 
mean-motion resonance perturbing planet, the libration 
period would be ~ 38.7 cycles and thus we do not possess 
sufficient temporal baseline to look for such perturbers 
yet. However, a 3:2 resonance would produce a libra- 
tion period of 15.4 cycles and thus we would expect the 
effects to be visible over the 12 observed cycles so far. 
This therefore excludes a 3:2 resonant perturber of mass 

> O.79M0 to 3-(T confidence. 

The current TTV data rules out an exomoon at the 
maximum orbital radius (for retrograde revolution) of 

> IO.GMq. The current TDV data rules an exomoon of 

> 17.3M0 at a minimum o rbital distance. 

Using the expressions of iFord fc HoW] ([20071) and 
assuming Mp » MTrojan, the expected libration period 
of Kepler-5b induced by Trojans at L4/L5 is 10.0 cycles 
and therefore we can search for such TTVs. We find 
such Trojans of cumulative mass > 3.14M0 at angular 
displacement 10° are excluded by our timings. Inspec- 
tion of the folded photometry at ±P/6 from the transit 
centre excludes Trojans of effective radius > 1.13i?® to 
3-a confidence. 

5.2.2. Periodograms 

An F-test periodogram of the TTV reveals only one sig- 
nificant peaks which is very close to Nyquist frequency. 
We therefore disregard this frequency as being signifi- 
cant. 

The TDV F-test periodogram reveals a significant peak 
at (4.0 ± 0.2) cycles with amplitude (37 ± 11) seconds 
with significance 2.2-ct. We again note that the period 
is not in the range which would be expected due to the 
known physical sources. The phasing effect has a period 
of 3 cycles and thus a mixing between the phasing and 
Nyquist frequency again offers a reasonable explanation 
for this signal, especially since peaks occur at 2 and 3 
cycles in the TDV periodogram as well. 

5.3. Depth and OOT Variations 

The transit depth is consistent with the globally fitted 
value yielding a of 7.9 for 11 degrees of freedom. Sim- 
ilarly, the OOT levels are flat yielding a of 4.0 for 11 
degrees of freedom. 

6. KEPLER-6b 

6.1. Global Fits 

6.1.1. Comparison to the discovery paper 

Kepler-6b was discovered bv lDunham et ahl ()2010[ ). In 
a similar manner to the other planets, we find consis- 
tent values for most parameters when compared to the 
iDunham et al.l (|2010[ ) analysis. The fitted models on the 



phase-folded lightcurves for method B are shown on Fig- 
ure [12] Correlated noise was checked for in the residuals 
using the Durbin- Watson statistic, which finds d = 1.943 
which is well within the 1% critical level of 1.870. The 
orbital fits to the RV points using method B are shown 
in Fig. I13[ depicting both the circular and eccentric fits. 
The final planetary parameters are summarized at the 
bottom of Table [61 

6.1.2. Eccentricity 

The global circular fit yields a — 9-01 and the ec- 
centric fit obtains = 7.75 for the RVs. Based on an 
F-test, the inclusion of two new degrees of freedom to 
describe the eccentricity is justified with a significance of 
< 1-(T The system is therefore best described with a cir- 
cular orbit with the current observations. We constrain 
e < 0.13 to 95% confidence. 

6.1.3. Secondary eclipse 

We here only consider the circular fit since the eccen- 
tric solution is not statistically significant. We find no 
evidence for a secondary eclipse detection for Kepler-6b. 
However, our analysis does exclude secondary eclipses of 
depth 51.5 ppm or greater to 3-(T confidence, which ex- 
cludes a geometric albedo > 0.32 to the same level. This 
therefore excludes the possibility of a planet covered in 
reflective clouds. A brightness temperature of 2445 K is 
ruled out to the same level, which places no constraints 
on the redistribution. 

6.1.4. Properties of the parent star Kepler-6 

The Yonsei-Yale model isochrones from IYi et al.l ([2001[ ) 
for metalhcity [Fe/H]=-|-0.34 are plotted in Figure [HI 
with the final choice of effective temperature Teff,^ and 
a/ Ri, marked, and encircled by the la and 2a confidence 
ellipsoids, both for the circular and the eccentric orbital 
solution. Here the MCMC distribution of a/i?* comes 
from the global modeling of the data. 

6.2. Transit timing analysis 
6.2.1. Analysis of variance 

We find TTV residuals of 17.4 seconds and TDV resid- 
uals of 25.4 seconds. The TTV indicates timings consis- 
tent with a linear ephemeris, producing a x^ of 7.8 for 11 
degrees of freedom. The TDV is also consistent with a 
non- variable system exhibiting a of 11.4 for 12 degrees 
of freedom. 

For the TTV, we note that for 11 degrees of freedom, 
excess scatter producing a x^ of 28.5 would be detected 
to 3-cr confidence. This therefore excludes a short-period 
signal of r.m.s. amplitude > 26.2 seconds. Similarly, for 
the TDV, we exclude scatter producing a x^ of 30.1, or 
a short-period r.m.s. amplitude of > 31.5 seconds, to 3- 
a confidence. This constitutes a 0.53% variation in the 
transit duration. 

These limits place constraints on the presence of per- 
turbing planets, moons and Trojans. For the case a 2:1 
mean-motion resonance perturbing planet, the libration 
period would be ~ 76.5 cycles and thus we do not possess 
sufficient baseline to look for such perturbers yet. How- 
ever, a 4:3 resonance would produce a libration period 
of 17.7 cycles and thus we would expect the effects to be 
visible over the 14 observed cycles so far. This therefore 
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Fig. 11. — Upper Left: TTV (squares) and TDV (diamonds) for Kepler-5b (see ^2.3.1 for details). Upper Right: TTV periodogram 
(solid) and TDV periodogram (dashed) for Kepler-5b, calculated using a sequential F-test (see ^^2.3.2 for description). Lower Left: 
Transit depths from individual fits of Kepler-5b (see %2.3.1 for details). Lower Right: "Phasing" of Kepler long cadence photometry for 
Kepler-5b (see ^2.3.4 for description). 



TABLE 5 

Mid-transit times, transit durations, transit depths and out-of-transit (normalized) fluxes for Kepler-5b. 



Epoch tc/(BJD-2,454,000) 



T/s 



Foot 



+U.UUU38 
0.00039 
+0.00039 
0.00040 
+0.00039 
0.00040 
+0.00039 
0.00040 
+0.00039 
-0.00039 
+0.00038 
0.00038 
+0.00041 
0.00041 
+0.00037 
0.00037 
+0.00038 
0.00038 
+0.00041 
0.00040 
+0.00039 
0.00038 
+0.00038 
0.00038 



D + '/y.U 



6S2+" "^^ 

U.DOZ_Q QQg 
U.OCiZ_Q QJ^J^ 

n fiS2+°"2'' 

U.OCiZ_Q (jj^j^ 
n KCJ7+0.022 

n fj57+0.033 

U.OZD_Q QQg 

U.D^O_Q QQg 
(,07 + 0.027 

U.DCil_Q QJ^Q 

U.DCiD_Q (jgg 



0.999952+^ 
0.999958t° 
0.999943t° 

'^0 
0.999970t° 

1+0 
.(+0 



UUUUi5 
000015 
000015 
000015 
000019 
000019 
000015 
000015 

.000015 
000015 
000015 
000015 
000015 
000015 
000015 
000015 

.000015 
000015 
000015 
000015 

.000015 
000015 
000015 
000015 




1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 



955.90106 
959.45010 
962.99799 
966.54678 
970.09481 
973.64409 
977.19164 
980.74048 
984.28922 
987.83677 
991.38619 
994.93424 



15177. 



15351.41^2.8 
15372.lt^«;i 



s+72.9 
-74.0 
3+73.6 
-71.9 
^+73. 5 
-72.5 

15284.6tg7g 

s+70:6 
-69.5 



15267.6: 
15241.91 
15222.4"* 



15176.5 
15231. 7t7i;i 
15242 
15392.9 
15239.6 



+77.7 
79.0 
+73.7 

72.9 
+71.5 
70.3 



0.999965"' 



0.999970"* 



0.999964_(, 

0.999966t° 

'^0 
0.999977t° 
-1+0 



0.999965"* 



0.999975"* 



0.999974"* 



excludes a 4:3 resonant perturber of mass > 0.38Af^ to 
3-cr confidence. 

The current TTV data rules out an exonioon at the 
maximum orbital radius (for retrograde revolution) of 

> 4.8M0. The current TDV data rules an exomoon of 

> 8.2M0 at a minimum or bital distance. 

Using the expressions of iFord fc HoW] (l2007h and 
assuming Mp >• Mxrojan, the expected libration period 
of Kepler-6b induced by Trojans at L4/L5 is 16.7 cycles 



and therefore we can search for such TTVs. We find 
such Trojans of cumulative mass > Q.QIM,^ at angular 
displacement 10° are excluded by our timings. Inspec- 
tion of the folded photometry at ±P/6 from the transit 
centre excludes Trojans of effective radius > O.STi?® to 
3-cr confidence. 



6.2.2. Periodograms 
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TABLE 6 

Global fits for Kepler-6b. Quoted values are medians of MCMC trials with errors given by 1-a quantiles. Two independent methods are 
used to fit the data (A&B) with both circular and eccentric modes (c&e).* = fixed parameter; f = parameter was floated but not fitted. 



Parameter 



Discovery 



Method Ax 



Method A.e 



Model B.c 



Model B.e 



Fitted params. 
P/days 
E (BJD-2,454,000) 

7^1,4/8 

T/s 

(Tl,2 ~ T3,4)/S 

(C/i?*)/days-i 
(Fp,d/F*)/ppm 

esinoj 

ecoso) 

7/ms~-' 

7/ms~-'-day~-'- 
itroj /days 
B 



q 0Q/17O1 +(I.(IU()U43 

954.48639lg;roil 
129821^1° 

1172011*2 

+234 
-84 
+0.047 
0.028 
+0.137 



O Qo^yr,, +0.000028 

954.48637lg:ri 
129561^^^ 
1172811°° 
I2OOI7P 

r,+0.036 
-0.021 



3.234723 ± 0.000017 
954.48636 ± 0.00014 



0* 
0* 

80.9 ± 2.6 
-18.3 ±3.5 
0* 
0* 

1.033 ± 0.004 



1211 
0.913 



0.069l„„e2 
14.801";- 

r;+36 
^-36 

0* 

0* 

78.4l«:J 
-18.2ltt 
0* 

-n 04i+°-°*2 

U.U41_Q 

1.033 ± 0.004 t 



0.910: 
0.061 



+0.107 
-0.054 

14.80l0;1^3 

-4l^2l 
042+° °^^ 

-0.0191°;°?^ 

-^'■'-3.0 
0* 

0* 

1.033 ± 0.004 t 



3.234720 ± 0.000018 
954.48640 ± 0.00014 
13003 ± 60 
11871 ± 16 
1227 ± 69 
0.901 ±0.011 
0.088l°:°« 
14.67 ± 0.02 
15± 18 
0* 
0* 

78.2 ±3.7 
-15.9 ±2.9 
0* 
0* 

1.033 ± 0.004 t 



3.234721 ± 0.000018 
954.48641 ± 0.00015 
12977 ± 60 
11762 ± 18 
1218 ±69 
0.901 ±0.011 

U.UOO_Q Q44 

14.69 ± 0.02 

21 ±35 
0.117 ±0.036 
012+° °^^ 
86.6 ± 1.6 
-13.6 ±2.7 
0* 
0* 

1.033 ± 0.004 t 



Tefl/K 


5647 ± 44 


5647 ± 88 


5647 ± 88 


5647 ± 88 


5647 ± 88 


log(g/cgs) 


4.59 ±0.10 


4.59 ±0.10 


4.59 ±0.10 


4.59 ±0.10 


4.59 ±0.10 


[Fe/H] (dex) 


+0.34 ± 0.08 


+0.34 ± 0.08 


+0.34 ± 0.08 


+0.34 ± 0.08 


+0.34 ± 0.08 


usini (kms~l) 


3.0 ±1.0 


3.0 ± 1.0 


3.0 ±1.0 


3.0 ±1.0 


3.0 ± 1.0 


Ml 




0.551-.1? 


0.54410:^2 


0.4471* 


0.4471* 


«2 




0.0ll°;f, 


o.o2i°:i^ 


0.2397* 


0.2397* 





0.0982910;^,^;^ 


nq^^+" ""^* 
u.uyoo_Q QQj^g 


0954+" ""l** 


0.0257 ± 0.0009 


0.0256 ± 0.0006 


a/Ri, 


7.05lg;JJ 


r, or, + 0.22 

'•■^^-0.48 


7 no+°-** 

(.UU_Q 4g 


7.21 ±0.18 


6.42 ± 0.26 


b 


0.398l«:g^^ 


26+"- 1" 


0.25l»-« 


0.296l»;»«« 


"•^"^-0.102 


ir 


86.8 ± 0.3 


87.911;* 


87.911;* 


87.7l°:« 


87.0 ±0.9 


e 


0* 


0* 


056+° °** 

U.UOD_Q 


0* 


0.120 ±0.036 


u;/° 






iir;+86 
J.J^O_24 




83 ±12 


RV jitter (ms-i) 




0.6 


0.0 


3.2 


2.7 




0.634 


o.709i;5:?«I 


0.6i9l»;lS 


0.68 ± 0.05 


0.48 ± 0.06 


log(9p/cgs) 


9 074+0.016 


3.0l0l0;°|? 


2.99l|;°|* 


3.01 ±0.03 


2.95 ±0.04 


Si,4/s 




12982+llg 


15371_2332 


13003 ± 60 


10022 ± 829 


tsec/(BJD-2,454,000) 




959.33827iro°o°oll 


957.68lJ;«° 


978.7470 ± 0.0001 


978.772 ± 0.048 



Derived stellar params. 



M^/Mq 


1 209+"-"** 


-0.065 


1 iqy+O.UbV 

-0.076 


1 Tif;+U.U34 

1.10D_Q Qy2 


^ ^-Q+0.064 
-^•-^'^-0.087 


R*/R@ 


1 0Q1 +0.017 
1-^91-0.034 


l.C)UD_g Q4g 


1 •370+0.109 

i.OIU_g ggj^ 


1.325 ±0.042 


1.503 ± 0.073 


log(9/cgs) 


4.236 ±0.011 




A 91 0+0. 050 
4-^l«-0.057 


4.24 ±0.02 


4.15 ±0.03 




j..a»_g 21 


, t:7+0.26 
-^•'^'-0.17 


1 71+0.32 
-^■'-^-0.25 


1.61 ±0.17 


2.06 ±0.28 


Mv 1 mag 




4 Q5+0.14 


^•^"-0.19 


4.33 ±0.13 


4.05 ±0.16 


Age/Gyr 


3.8 ±1.0 




5 65+2.84 




r; 7+2.6 
'^■'-1.0 


Distancc/pc 




604is° 


630l5i 


621 ± 22 


704 ± 37 



Derived planetary params. 



0.025 
-0.030 



+0.U52 
0.051 



645+"-"*''' 
U.D+0_Q Q44 

1 971+0.116 

^-0.091 

0.389l°;°|t 



Mp/Mj 
RpIRi 

a/AU 



0.669 
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1 323+°° 

1.0Z0_Q Q29 

n QC9+0.018 

U.OOZ_g Q22 

04^67+°-°°°'^'^ 
1500 ± 200 



1 9ns+oi29 

l.^U»_Q 049 

n 426+0 069 
^■^^"-O.IOB 



0.04452l°;roi 



1480"* 



0.04466l°;°°°51 

'til 



1511"' 



0.616 ±0.035 
1.224 ±0.043 
0.41 ± 0.04 

0447+° °°°* 
1488 ± 30 



0.692 ±0.032 
1.389 ±0.072 
0.32 ± 0.04 

0451+° °°°* 

U.U401_o 0012 

1579 ± 41 
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Fig, 12, — Top: Phase-folded primary transit lightcurve of 
Kepler-6for method B, The upper curve shows the circular fit, the 
bottom curve the eccentric fit. Solid (blue) lines indicate the best 
fit resampled model (with bin- number 4). The dashed (red) lines 
show the corresponding unbinned model, which one would get if the 
transit was observed with infinitely fine time resolution. Residuals 
from the best fit resampled model for the circular and eccentric 
solutions are shown below in respective order.. 
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Phase with respect to T^, 

Fig. 13. — Top: RV measurements from Keck for Kepler-6, 
along with an orbital fit, shown as a function of orbital phase, us- 
ing our best-fit period. Solid (blue) line shows the circular orbital 
fit with binned RV model (3 bins, separated by 600 seconds). The 
dashed (red) line shows the eccentric orbital fit (with the same 
bin parameters). Zero phase is defined by the transit midpoint. 
The center-of-mass velocity has been subtracted. Note that the er- 
ror bars include the stellar jitter (taken for the circular solution), 
added in quadrature to the formal errors given by the spectrum 
reduction pipeline. Fits from method B. Middle: phased residu- 
als after subtracting the orbital fit for the circular solution. The 
r.m.s. variation of the residuals is S./ms^-"^, and the stellar jitter 
that was added in quadrature to the formal RV errors is 3.2 m s~^. 
Bottom: phased residuals after subtracting the orbital fit for the 
eccentric solution. Here the r.m.s. variation of the residuals is 
5.5ms~^, and the stellar jitter is 2.7ms~-'^. 

The TTV periodogram reveals a significant peak 
of period (5.34 ± 0.26) cycles with amplitude (19.7 ± 
5.0) seconds and FAP 2.6-a, which is one of our high- 
est significances found, yet still below our formal detec- 
tion threshold of i-a. Nevertheless, the signal warrants a 
more detailed investigation as it does not appear to land 
on any of the probable aliasing frequencies. 




11.0 I ' ' ' ' ' ' ' ' ' ' 1 

6000 5500 5000 
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Fig. 14. — Stellar isochrones from lYi et al.l 1 120011) for metal- 
licity [Fe/H]=+0.34 and ages 1.0, 1.4, 1.8, 2.2, 2.6, 3.0, 3.4 and 
3.8Gyrs. The final choice of Tcfj* and a/R* for the circular solu- 
tion are marked by a filled circle, and is encircled by the Icr and 
2(T confidence ellipsoids (solid, blue lines). Corresponding values 
and confidence ellipsoids for the eccentric solution are shown with 
a triangle and dashed (red) lines. Fits from method B. 

Given that Trojans should induce libration of pe- 
riod 16.7 cycles, this hypothesis can be instantly dis- 
counted. For two planets in a j:j + 1 mean motion- 
resonance, we may set the TTV signal's period to be 
equal to th e expressio n for t he expected libration pe- 
riod (from lAgol et all (|2005[ )) and solve for j. This 
yields j = 7.37l[];^^, or a period ratio of 1.136 ± 0.005, 
which does not appear to be a resonant configuration. 
IBarnes k Greenbergj (|2006l ) used the Hill criterion to 
provide an expression for estimating period ratios which 
are not guaranteed to be dynamically stable (equation 
2). Using this expression, we estimate period ratios of 
> 0.75 and < 1.34 are unlikely to be stable for masses 
above IM^. This makes the hypothesis of a mean- motion 
resonance unlikely. 

For an out-of-resonant signal, the period of the TTV 
is more variable, depending on the difference between 
the perturber's period and the nea rest reson ant posi- 
tion, quantified by the parameter e in lAgol et" al. (2005). 
For a completely non-resonant system, e = 1 and the 
period ratio is 6.34. As we move closer to resonance, 
this period ratio decreases. For the maximally non- 
resonant configuration, the perturber would require a 
mass of Afport — (0.75 ± 0.19)Mj to generate the pu- 
tative TTV signal. Given that this is approximately the 
same mass as the transiting planet and the period ratio is 
not large, it is clear that such a signal would be evident 
in the RV data. This makes the hypothesis of a non- 
resonant perturber unlikely but not completely excluded 
due to the possibility of ou t -of-the-pla ne perturbers (see 
iNesvorny fc Bcaugd IMS); l^raTeTal. (2010)). 

If such a signal was due to an aliasing of an exomoon's 
orbital period, the maximum amplitude would occur for 
a moon at the maximum orbital separation given by 
iDomingos et al.l ()2006[) . for a retrograde configuration. 
A moon of mass (1.8±O.4)M0 could generate such a sig- 
nal and closer moons require larger masses. Such a moon 
should produce a photometric dip of between 40ppm and 
170ppm, depending on its composit ion. Following the 
method outlined in iKippind (|2009cO . we may use the 
phasing predicted by the TTV fit to folded the individual 
transit residuals for cases where maximum planet-moon 
separation is predicted; epochs 1, 7 & 11 and the time- 
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reversed epochs 4 & 9. Inspection of this folded residual 
rules out exomoons of i? > 1.6i?© to 2-a confidence. 

An additional method of confirming a moon would 
be through a TDV analysis. The TDV periodogram 
does have one significant peak at (3.6 ± 0.2) days though 
but with a weaker significance of 2. 2- a and amplitude 
(24.3 ± 7.7) seco nds. Using the ratio of the TTV and 
TDV amplitude, iKippinsj (|2009al) showed a solution for 
the exomoon period is possible. Following this method 
and assuming both the TTV peak and the TDV peak are 
genuinely due to a moon, we estimate an exomoon period 
of (4.2± 1.7) hours which would place the moon at a semi- 
major axis of (0.8 ± 0.2) planetary radii, i.e. within the 
planet's volume. Although there is a large uncertainty on 
this, we conclude that the hypothesis of an exomoon in- 
ducing both the TTV and TDV putative signals is highly 
improbable. 

The TDV peak is quite likely to be spurious due to the 
proximity of the periodogram peak to twice the Nyquist 
frequency. The exomoon hypothesis can be further in- 
terrogated by evaluating the range of dynamically sta- 
ble moons around Kepler-6b. Assuming Jovian-like tidal 
dissipation values of fc2p = . 5 and Qp — 10^, the ex- 
pression of lBarnes fc O'BrienI ()2002D indicate that a sin- 
gle moon could only be stable for the age of the system 
if it has a mass below O.25M0. In contrast, the TTV 
could only be produced by moons > l.SM^ suggesting 
we would require Qp > 10^ to make the two compati- 
ble. In conclusion, there is no confirming evidence for 
the moon hypothesis and at thus we cannot assign this 
explanation to the TTV peak. 

The TTV signal's period of 17.3 days is similar to that 
what we might expect for the star's rotational period and 
thus could be due to the motion of star spots across the 
stellar surface, as observed by Alonso et al, (2008) for 
CoRoT-2b. Based upon the projected vsmi rotation of 
the star and the inferred stellar radius, we estimate the 
rotation period to be 23.515^9^^ days which is compatible 
with the period of the TTV signal. A bisector analysis 
is not possible due to the very low cadence but will be 
possible once short-cadence data becomes available. We 
propose that this is the most likely explanation for the 
signal, based upon the current evidence. 

The phasing signal has a period of 4 cycles, which is 
close to both periods found in the TTV and TDV pe- 
riodograms. This therefore detracts credence to the hy- 
pothesis that these signals are genuine. Further observa- 
tions, particularly in short-cadence mode, where phasing 
will not be a significant issue, will clarify the reality of 
these signals. 

6.3. Depth and OOT Variations 

The transit depth is consistent with the globally fitted 
value yielding a of 10.1 for 12 degrees of freedom. 
Similarly, the OOT levels are flat yielding a of 5.1 for 
12 degrees of freedom. 

7. KEPLER-7b 
7.1. Global fits 
7.1.1. Comparison to the discovery paper 

Kepler- 7b was discovered by iLatham et al.l ()2010[ ). 

The planet is noteworthy for its very low density and 
thus significant inflation. Our global fits find highly con- 



sistent r esults with the origin al discovery paper. We note 
that the ILatham et all ()2010|) values are the most consis- 
tent with our own values out of any of discovery papers 
discussed in this work. The fitted models on the phase- 
folded lightcurves are shown on Figure [TBI for method B. 
Correlated noise was checked for in the residuals using 
the Durbin- Watson statistic, which finds d ~ 1.866, at 
the border of the 1% critical level of 1.868, suggesting 
correlated noise is marginal. The orbital fits to the RV 
points are shown in Fig. I19[ depicting both the circu- 
lar and eccentric fits using method B. Unusually, both 
methods A and B found stellar jitter levels consistent 
with zero, which suggests the publicly available RV mea- 
surements have perhaps overestimated errors. The final 
planetary parameters are summarized at the bottom of 
Table [51 

7.1.2. Eccentricity 

The global circular fit of method A yields a = 6.87 
and the eccentric fit obtains = 4.89 for the RVs. An 
F-test finds that the eccentric fit has a significance of 
79.4%. Therefore, we conclude that the eccentric fit does 
not meet our detection criteria and we constrain e < 0.31 
to 95% confidence. 

One obvious question at this point is whether the range 
of allowed eccentricities could explain the bloated radius 
of Kepler- 7b through tidal heating. The Qp value of the 
planet is an element of large uncertainity here but we can 
constrain Q's minimum limit using the same method we 



+ 14.0 
-4.6 



adopted for Kepler-4b. We find that Qp > 9.4 
10^. One may argue that this is only valid under the 
assumption that the eccentric model is selected over the 
circular one, but we would counter that in the circular 
model case the tidal heating is zero anyway and thus the 
hypothesis of tidal heating being the origin of the bloated 
radius is invalid. 

Using the posterior distribution of the minimum al- 
lowed Qp value, we can translate this to a maximum 
allowed tidal en e rgy d issipation using the expressions of 
IPeale fc CassCTl (fl978l) : 
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GM^m 



2 Q 



P.min 



(16) 



We divide this value by the amount of power hit- 
ting the planetary surface, found by L*i?p/(4a^) to give 
us the posterior of the fraction of maximum allowed 
tidal heating relative to the incoming radiation and find 
^^tidai/^'stciiar — ^-^tf o X 10~^^. We therefore conclude 
that tidal heating plays a negligible role in the heating 
of Kepler- 7b and the planet's inflated radius. 

7.1.3. Secondary eclipse 

As before, method B is the better tool to constrain 
secondaries due to its rescaling method. Method B de- 
tects a secondary eclipse in the circular fits of depth 
(47 l b 14)ppm and s ignific ance 3.5-a (99.96%). We note 
that ILatham et al.l (|2010l ) reported a 2.4-cr weakly de- 
tected eclipse but did not report a depth estimate. The 
MCMC distribution of the Fp^d/Fi, values for method B 
is shown in Figure [T51 

Inspection of the method A and ILatham et al.l ()2010f ) 
values relative to the method B result reveals a clear dif- 
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Fig. 15. — Upper Left: TTV (squares) and TDV (diamonds) for Kepler-6b (see ^2.3.1 for details). Upper Right: TTV periodogram 
(solid) and TDV periodogram (dashed) for Kepler-6b, calculated using a sequential F-test (see ^2.3.2 for description). Lower Left: 
Transit depths from individual fits of Kepler-6b (see %2.3.1 for details). Lower Right: "Phasing" of Kepler long cadence photometry for 
Kepler-6b (see %2.3.4 for description). 



TABLE 7 

Mid-transit times, transit durations, transit depths and out-of-transit (normalized) fluxes for Kepler-6b. 



Epoch tc/(BJD-2,454,000) 
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11758. Slgjt 

h50!8 
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999971 +0 000015 

-,+0.000015 
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ference and we must consider why we obtain such differ- 
ent numbers. All of the fits are circular and thus the or- 
bital solutions are nearly identical with the eclipse occur- 
ing at the same phase in each case. The only difference 
between them is the definition of the out-of-eclipse base- 
line. Method A assumes the same baseline level for both 
primary and secondary transit whereas method B uses a 
simplified phase curve model as described in §2.2. If the 
secondary eclipse local baseline is different from the pri- 
mary transit baseline, method A would move towards an 
average solution and actually attenuate the eclipse depth. 
The s ame attenuation also occurs in the ILatham et al.l 
(|2010D which utilizes the same technique as method A 
(D. Latham, personal communication). 

To investigate this futher, we may compare the three 
baseline levels found by method B: Foot = (0.999976 ± 
0.000007), Foorr, = (0.999972 ± 0.000003) and Foos = 
(0.999993± 0.000007). The first two values are compati- 
ble with a constant but the baseline surrounding the sec- 
ondary eclipse is clearly larger. If the secondary eclipse 
baseline is larger, method A would favor a lower average 
baseline due to Foot and Foom and thus find a lower sec- 
ondary eclipse. This seems to be the explanation with 
method A finding Foot = (0.9999764 ± 0.0000024). We 
believe this explains the discrepancy between the two 
methods. 

Using the method B results only, the difference in base- 
lines indicates a day-night contrast. Using our MCMC 
trials, we estimate a nightside fiux (Fn/i^*) — 29^11 ppm 
which is of 2.4-cr significance. We can also calculate the 
day- night difference explicitly and we find [(Fp^d/-F*) — 
{Fp,n/F*)] = (17 ± 9) ppm. Since nightside flux cannot 
be due to reflected light, this supports the hypothesis 
of thermal emission as a source for both the secondary 
eclipse and the phase curve. 

There is a very strong caveat to bear in mind for the 
day-night contrast. The data we have analyzed has been 
processed by the Kepler pipeline which involves numer- 
ous data manipulation processes. Time trends of very 
low amplitude and of time periods of days may be af- 
fected by this processing, with the most probable be- 
haviour being that of a long-cut fllter. As a result, a 
strong-day night contrast could be attenuated by the act 
of passing through this filter. It is therefore still possi- 
ble that the day-night contrast is 100% (i.e. equal to the 
secondary eclipse depth) and thus the light we are see- 
ing is reflected light, rather than thermal emission, but 
the phase curve is artiflcially distorted by the pipeline 
processing. 

If the secondary eclipse was due to thermal emission, 
we may estimate the brightness temperature by integrat- 
ing the Planck function of the planet over the Kepler re- 
sponse function, divided by that of the star. We repeat 
this for planetary brightness temperatures from 2000K 
to 3000K in IK steps and calculate (Fpd/-F*) in each 
case. We find that a best-fit brightness temperature of 
2570j^g5^K, which is much larger than the equilibrium 
temperature of (1554 ± 32) K. Even the nightside flux 
of (Fp^n/Fi,) = 29^1^4 ppm would require a temperature 

for the nightside of 2441^i[7gK. These very large tem- 
peratures do not seem consistent with the energy bud- 
get of the planet, particularly as tidal heating has al- 
ready been shown to offer a negligible contribution to 



the energy budget of this planet. From stellar irradia- 
tion only, the redistribution factor would have to be 7.48 
to re-create the dayside eclipse, which is greater than 
the maximum physically permitted value of 8/3. The 
only remaining possibility is internal heating, for which 
we find a heat source responsible for 2440 K could ex- 
plain both the dayside and nightside fluxes. This would 
require an internal heat source of \og{L/LQ) = —3.19, 
which is equivalent to the luminosity of an M6-M7 star. 
For example, LHS 292 (M6.5) and VB8 (M7) have lumi- 
nosities of log(L/io) = —3.19 and log(i/L0) = —3.24 
respectively. 

Whilst other sources of heating are possible apart from 
tidal, for example radioactivity, the magnitude required 
detracts credence from the hypothesis of thermal emis- 
sion as the principal origin for the eclipse depth. 

If due to reflection, method B flnds a geometric albedo 
of Ag = (0.38 ± 0.12), which is somewhat larger than 
the value found for Kepler-5b of Ag = (0.15 ± 0.10), but 
not greatly different from the albedos of Jupiter and Sat- 
urn (0.52 and 0.47 respe c tively ) and consistent with the 
study of iCowan fc AgoU ()2010l ). We believe this latter 
hypothesis, as with Kepler-5b, is a more plausible sce- 
nario for the secondary eclipse measurement. As a con- 
sequence, we should not expect any nightside fiux, which 
bears some questions on the apparent non-zero nightside 
flux detected in this study. There are two possible ex- 
planations which are consistent with reflection: i) the 
nightside flux is not statistically significant at 2.4-(t ii) 
the phase curve has been attenuated due to the Kepler 
pipeline producing effects which mimic a long-cut filter 
iii) Foos is increased away from zero nightside flux due 
to the phase curve of the planet i.e. refiected light from 
the viewable crescent. In reality, a combination of these 
seems probable. Once more transits become available for 
analysis, a more sophisticated phase curve model can be 
fitted to the data set to resolve this question. 
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Fig. 16. — Top: Phase-folded primary transit lightcurve of 
Kepler-7 using method B. The upper curve shows the circular fit, 
the bottom curve the eccentric fit. Solid (blue) lines indicate the 
best fit resampled model (with bin- number 4). The dashed (red) 
lines show the corresponding unbinned model, which one would 
get if the transit was observed with infinitely fine time resolution. 
Residuals from the best fit resampled model for the circular and 
eccentric solutions are shown below in respective order. 



7.1.4. Properties of the parent star Kepler-7 
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Fig. 17. — Phase-folded secondary transit lightcurve of Kepler-7 
using method B (top), and residuals from the best fit (bottom). 
Only the fit for the circular orbital solution is shown. 
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Fig. 18. — Distribution of Fp^d/F* for for Kepler-7 from the 
global modeling of method B. Fp ,i/F-i, is primarily constrained by 
the depth of the secondary eclipse. Only results for the circular 
orbital solution are shown. 

The Yonsei-Yale model isochrones from lYi et al.l (|2001[ ) 
for metallicity [Fe/H]=-|-0.11 are plotted in Figure 
with the final choice of effective temperature Tcff^, and 
a/ Ri, marked, and encircled by the Icr and 2fj confidence 
ellipsoids, both for the circular and the eccentric orbital 
solution. Here, the MCMC distribution of a/R^, comes 
from the global modeling of the data. 

7.2. Transit timing analysis 
7.2.1. Analysis of variance 

We find TTV residuals of 32.1 seconds and TDV resid- 
uals of 22.2 seconds. The TTV indicates timings con- 
sistent with a linear ephemeris, producing a of 9.1 
for 7 degrees of freedom (75.1% significance of excess). 
The TDV is also consistent with a non-variable system 
exhibiting a of 4.1 for 8 degrees of freedom (84.9% 
significance of being too low). 

For the TTV, we note that for 7 degrees of freedom, 
excess scatter producing a of 21.8 would be detected 
to 3-cr confidence. This therefore excludes a short-period 
signal of r.m.s. amplitude > 38.0 seconds. Similarly, for 
the TDV, we exclude scatter producing a x^ of 23.6, or 
a short-period r.m.s. amplitude of > 42.6 seconds, to 3- 
cr confidence. This consitutes a 0.50% variation in the 
transit duration. 

These limits place constraints on the presence of per- 
turbing planets, moons and Trojans. For the case a 2:1 
mean-motion resonance perturbing planet, the libration 
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Fig. 19. — Top: I{V measurements from Keck for Kepler-7, 
along with an orbital fit, shown as a function of orbital phase, us- 
ing our best-fit period. Solid (blue) line shows the circular orbital 
fit with binned RV model (3 bins, separated by 600 seconds). The 
dashed (red) line shows the eccentric orbital fit (with the same 
bin parameters). Zero phase is defined by the transit midpoint. 
The center-of-mass velocity has been subtracted. Note that the er- 
ror bars include the stellar jitter (taken for the circular solution), 
added in quadrature to the formal errors given by the spectrum 
reduction pipeline. Fits from method B. Middle: phased residu- 
als after subtracting the orbital fit for the circular solution. The 
r.m.s. variation of the residuals is 8.1 ms~^, and the stellar jitter 
that was added in quadrature to the formal RV errors is 0.0 m s~^. 
Bottom: phased residuals after subtracting the orbital fit for the 
eccentric solution. Here the r.m.s. variation of the residuals is 
6.9 ms^^, and the stellar jitter is 0.0 ms~^. 




Effective temperature [K] 

Fig. 20. — Stellar isochrones from lYi et al.l 1 120011) for metal- 
licity [Fe/H]=+0.11 and ages 1.0, 1.4, 1.8, 2.2, 2.6, 3.0, 3.4 and 
3.8Gyrs. The final choice of T^gi, and afRt, for the circular solu- 
tion are marked by a filled circle, and is encircled by the Icr and 
2(T confidence ellipsoids (solid, blue lines). Corresponding values 
and confidence ellipsoids for the eccentric solution are shown with 
a triangle and dashed (red) lines. Fits from method B. 

period would be ~ 111 cycles and thus we do not possess 
sufficient baseline to look for such perturbers yet. After 
a further 102 cycles have been obtained at the same pho- 
tometric quality, the presence of a 2:1 resonant perturber 
of mass O.O6M0 may be constrained. 

The current TTV data rules out an exomoon at the 
maximum orbital radius (for retrograde revolution) of 

> 2.5M0. The current TDV data rules an exomoon of 

> 5.9M0 at a minimum orbital distance. 

Using the expressions of ,Ford fc HolmanI (|2007l) and 
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TABLE 8 

Global fits for Kepler-lb. Quoted values are medians of MCMC trials with errors given by 1-a 
used to fit the data (A&B) with both circular and eccentric modes (c&e).* = fixed parameter; \ 



quantiles. Two independent methods are 
= parameter was floated but not fitted. 



Parameter 



Discovery 



Method A.c 



Method A.e 



Model B.c 



Model B.e 



Fitted params. 
P/days 
E (BJD-2,454,000) 

7^1,4/8 

T/s 

(Tl,2 ~ T3,4)/S 

(C/i?*)/days-i 
{Fp,d/^*)/ppm 
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0.000061 
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0* 
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0.00 
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SME derived params. 



Teff/K 


5933 ± 44 
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5933 + 88 
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log(g/cgs) 
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Derived stellar params. 
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Derived planetary params. 
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a/AU 

TeJK 



0.42 
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0.1671°;°^^ 

0.0608l°;°gl^ 
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assuming Mp ^ AfTrojan, the expected libration period 
of Kepler-7b induced by Trojans at LA/Lb is 21.4 cycles 
and therefore we do not yet possess sufficient baseline to 
search for these TTVs. Inspection of the folded photom- 
etry at ±P/6 from the transit centre excludes Trojans of 
effective radius > l.lli?® to S-cr confidence. 

What is peculiar about this system is that the TTV 
produces a slight excess variance whilst the TDV does 
the opposite. The error on the TDV data points, and 
thus the scatter, should be exactly the same as the TTV 
points, or in fact slightly larger due to limb darkening 
effects. Thus, the observation of the TTV scatter being 
much larger than the TDV scatter is curious. We note 
for all other planets we studied we found the TDV scat- 
ter is greater than or equal to the TTV scatter. Also, 
for several planets we find the observed scatter is much 
lower than the uncertainties predict, suggesting we have 
overestimated our errors. If this is true for the Kepler- 
7b system, the TDV scatter indicates that the measure- 
ment uncertainties should be scaled by a factor of 0.71. 
This factor should be common to the TTV, which would 
change the observed from 9.1 to 17.7, for 7 degrees 
of freedom. The scaled would have a significance of 
98.7%, or 2.5-a, which is still below our formal detection 
threshold. Nevertheless, we identify Kepler- 7b as an im- 
portant target for continued transit timing analysis. 

7.2.2. Periodograms 

The TTV F-test periodogram presents one significant 
peak very close to Nyquist frequency, which we may 
discard. The TDV periodogram presents no significant 
peaks. 

7.3. Depth and OOT Variations 

The transit depth is consistent with the globally fitted 
value yielding a of 5.8 for 8 degrees of freedom. Sim- 
ilarly, the OOT levels are flat yielding a x^ of 3.4 for 8 
degrees of freedom. 

8. KEPLER-8b 

8.1. Global Fits 

8.1.1. Comparison to the discovery paper 

Kepler-8b was discovered by iJenkins et al.l ()2010aD . 
Kepler-8b posed several unique problems not encoun- 
tered with the other Kepler planets. The data format of 
the public data was quite different from the other sets. 
Typical values of the photometric data was ~ whereas 
the other data sets where ~ 1. This suggests that either 
the photometry was raw differential magnitudes or nor- 
malized fluxes with unity subtracted. We tried both as- 
sumptions and that the normalized fluxes minus unity as- 
sumption led to system parameters more consistent with 
the original discovery paper. This assumption has been 
since confirmed (personal communication with J. Jenk- 
ins). 

Another peculiarity of the data is that the r.m.s. scat- 
ter is not constant but actually increases significantly 
with time, as discussed in §3.2. The source of this varying 
scatter is an instrumental effect, caused by oscillations 
in the photometry as a result of a heater on one of the 
two reaction w heel assemblies that m odulates the focus 
by ~ 1 /xm fsee IJenkins et ahl (|2010cf ) for more details). 
Since the discovery of this effect in the Ql photometry. 



the Kepler Team have subsequently trimmed the bound- 
ing temperatures for this heater to reduce the amplitude 
and heater cycle (J. Jenkins, personal communication). 
Therefore, future data sets will not suffer from this effect 
as accutely as the photometry analyzed here. Concor- 
dantly, the results obtained in our analysis may not be 
as reliable as for the other Kepler planets. 

In our analysis, we disregard radial velocity points 
occuring during the transit to avoid the Rossiter- 
McLaughlin effect. This is justified since our prinicpal 
goal here is to produce refined system parameters rather 
measure the spin-orbit alignment. 

We find the orbit is consistent with a circular orbit 
and derive a significantly lower plan etary mass than that 
reported by IJenkins et al.l ()2010af) . The fitted models 
on the phase- folded lightcurves are shown on Figure [531 
The orbital fits to the RV points are shown in Fig. IMl 
depicting both the circular and eccentric fits. The final 
planetary parameters are summarized at the bottom of 
Table [lOl 

The residuals of the lightcurve fits from method A (fit- 
ted limb darkening) were inspected using the Durbin- 
Watson statistic, which finds d = 1.344, indicating time 
correlated noise is certainly present. The origin of this is 
very likely the heater oscillations described earlier. An 
F-test periodogram finds the largest peak occuring at 
3.27 hours with amplitude 50.6 ppm, significance 4.25-(t. 
The rich forest of pcrdiodic lines make a correction un- 
feasible, as seen in Figure [531 

8.1.2. Eccentricity 

The global circular fit yields a x^ = 20.96 and the 
eccentric fit obtains x^ — 20.55. Based on an F-test, 
the inclusion of two new degrees of freedom to describe 
the eccentricity is justified at a significance of 13.1% and 
thus is not accepted. 

8.1.3. Secondary eclipse 

Using the circular fit, there is no evidence for a sec- 
ondary eclipse; we can exclude secondary eclipses of 
depth 101.1 ppm to 3-cr confidence, or equivalently ge- 
ometric albedos > 0.63. We can also exclude a bright- 
ness temperature > 2932 K to the same confidence level, 
which places no constraints on the redistribution of heat 
around the planet. 

8.1.4. Properties of the parent star Kepler-8 

The Yonsei-Yale model isochrones from iYi et al.l ()2001[ ) 
for metallicity [Fe/II]=— 0.055 are plotted in Figure 1551 
with the final choice of effective temperature Tcff* and 
a/ Ri, marked, and encircled by the la and 2a confidence 
ellipsoids, both for the circular and the eccentric orbital 
solution. Here the MCMC distribution of a/i?* comes 
from the global modeling of the data. 

8.2. Transit timing analysis 
8.2.1. Analysis of variance 

We find TTV residuals of 68.1 seconds and TDV resid- 
uals of 92.2 seconds. The TTV indicates timings consis- 
tent with a linear ephemeris, producing a x^ of 10.7 for 
11 degrees of freedom. The TDV is also consistent with a 
non-variable system exhibiting a x^ of 18.8 for 12 degrees 
of freedom (84.9% significance of excess). 
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Fig. 21. — Upper Left: TTV (squares) and TDV (diamonds) for Kepler-7b (see ^2.3.1 for details). Upper Right: TTV periodogram 
(solid) and TDV periodogram (dashed) for Kepler-7b, calculated using a sequential F-test (see ^^2.3.2 for description). Lower Left: 
Transit depths from individual fits of Kepler-7b (see %2.3.1 for details). Lower Right: "Phasing" of Kepler long cadence photometry for 
Kepler-7b (see ^2.3.4 for description). 



TABLE 9 

Mid-transit times, transit durations, transit depths and out-of-transit (normalized) fluxes for Kepler-7b. 



Epoch tc/(BJD-2,454,000) 



T/s 



Foot 



957.50508±1^;S 

r|+0. 00035 
^-0.00034 
D+0.00034 



+64.8 



U.D4i_Q gj^4 

U.OOO_Q 
U.DD4_Q Q23 
U.Om_Q (J26 

U.D40_Q gj^g 

U.D40_Q Q20 

0.664^0.025 
0.680^0.020 
0.640to.oj8 



UUUUiU 
000010 
000011 
000011 
000010 
000010 
000010 
000010 
.000010 
000010 
000011 
000011 
000010 
000010 
000010 
000010 
000010 
000010 



962.39000"* 



967. 27558_Q 00033 

0+0.00034 
-0.00033 
-+0.00032 
-0.00032 
-,+0.00035 
-0.00034 
j+0. 00035 
-0.00035 
.+0.00038 
-0.00038 

996.58887^0.00033 



972.16172: 
977.04645^ 
981.93250^! 
986.81768^ 
991.70414+ 



17006.2 



-65.6 
S+65.0 
'-65. 1 

17066.2tg2.3 

16963.6tg2| 

7+61.> 
-60.7 
-,+61.9 
-61.2 
.,+63.3 
-62.2 
3+66.6 



16993.! 



16978.2: 
16922.01 
16948.6"* 



16951.3" 



65.7 



16931. 2t^^f 



0.999970 
0.999967 
0.999975 
0.999974 
0.999973 
0.999976 
0.999980 
0.999976 
0.999957 



For the TTV, we note that for 11 degrees of freedom, 
excess scatter producing a of 28.5 would be detected 
to 3-cr confidence. This therefore excludes a short-period 
signal of r.m.s. amplitude > 58.8 seconds. Similarly, for 
the TDV, we exclude scatter producing a of 30.1, or 
a short-period r.m.s. amplitude of > 73.0 seconds, to 3- 
<T confidence. This consitutes a 1.46% variation in the 
transit duration. 

These limits place constraints on the presence of per- 
turbing planets, moons and Trojans. For a 2:1 mean- 
motion resonance, the libration period would be ^ 
83 cycles and thus we do not possess sufficient baseline 



yet to look for such perturbers. However, a 4:3 resonance 
system would yield a libration period of ~ 19 cycles and 
therefore we would expect to see such a signal with the 13 
cycles observed. This means a 4:3 perturber of > O.5OM0 
is ruled out to 3-cr confidence by the current TTV data 
for this object. After 83 cycles have been obtained at the 
same photometric quality, the presence of a 2:1 resonant 
perturber of mass > O.ltM^ may be constrained. 

The current TTV data rules out an exomoon at the 
maximum orbital radius (for retrograde revolution) of 

> 2.1Af0. The current TDV data rules an exomoon of 

> 2I.5M0 at a minimum orbital distance. 
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TABLE 10 

Fits for Kepler-8b. Quoted values are medians of MCMC trials with errors given by l-a quantiles. Two independent methods are used 
to fit the data (A&B) with both circular and eccentric modes (c&e).* = fixed parameter; f = parameter was floated but not fitted. 



Parameter Discovery Method A.c Method A.e Model B.c Model B.e 

Fitted params. 



P/days 3.522541^;^^^^^ 3.522261o;ooqJ3 3.522260to oooog5 3.522418 ± 0.000057 3.522416 ± 0.000057 
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0* 


0* 


0* 


0* 


0* 
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0* 
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0* 
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0* 
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0* 
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SME derived params. 


Teff/K 


6213 ± 150 


6213 ± 150 


6213 ± 150 


6213 ± 150 


6213 ± 150 


log(s/cgs) 


4.28 ±0.10 


4.28 ±0.10 


4.28 ±0.10 


4.28 ±0.10 
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[Fe/H] (dex) 
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-0.055 ± 0.06 
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10.5 ±0.7 
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M2 
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Derived stellar params. 
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Derived planetary params. 
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Fig. 22. — F-test periodogram of residuals from primary transit 
lightcurve of Kepler-8. The y-axis gives the false-alarm-probability 
subtracted from unity, in sigmas confidence. The highest peak 
occurs for a period of 3.27 hours with amplitude 50.6 ppm. 
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Fig. 23. — Top: Phase-folded primary transit lightcurve of 
Kepler-8. The upper curve shoves the circular fit, the bottom curve 
the eccentric fit. Solid (blue) lines indicate the best fit resampled 
model (with bin- number 4). The dashed (red) lines show the corre- 
sponding unbinned model, which one would get if the transit was 
observed with infinitely fine time resolution. Residuals from the 
best fit resampled model for the circular and eccentric solutions 
are shown below in respective order. 



Using the expressions of iFord fc HoW] ([20071) and 
assuming Mp ^ A^Trojan, the expected hbration period 
of Kepler-8b induced by Trojans at L4/L5 is 25.5 cycles 
and therefore we do not yet possess sufficient basehne to 
search for these TTVs. Inspection of the folded photom- 
etry at ±P/6 from the transit centre excludes Trojans of 
effective radius > 1.58i?® to 3-cr confidence. 

8.2.2. Periodograms 

The TTV periodogram reveals a peak at twice the 
Nyquist period, which we disregard. No other signifi- 
cant peaks are present. The TDV also has no significant 
peaks. The phasing signal appears to have a period of 2 
cycles. 




Phase with respect to T^, 

Fig. 24. — Top: I{V measurements from Keck for Kepler-8, 
along with an orbital fit, shown as a function of orbital phase, using 
our best-fit period. Solid (blue) line shows the circular orbital fit 
with binned RV model (3 bins, separated by 600 seconds). The 
dashed (red) line shows the eccentric orbital fit (with the same bin 
parameters). Zero phase is defined by the transit midpoint. The 
center-of-mass velocity has been subtracted. Note that the error 
bars include the stellar jitter (taken for the circular solution), added 
in quadrature to the formal errors given by the spectrum reduction 
pipeline. Middle: phased residuals after subtracting the orbital 
fit for the circular solution. The r.m.s. variation of the residuals 
is 45.7ms~^, and the stellar jitter that was added in quadrature 
to the formal RV errors is 31.5 ms~^. Bottom: phased residuals 
after subtracting the orbital fit for the eccentric solution. Here the 
r.m.s. variation of the residuals is 42.8ms~'^, and the stellar jitter 
is 27.1 ms^-'^. 
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Fig. 25. — Stellar isochrones from lYi et al.l 1 120011) for metal- 
licity [Fe/H]=-0.055 and ages 1.0, 1.4, 1.8, 2.2, 2.6, 3.0, 3.4 and 
3.8Gyrs. The final choice of T^ffi, and a/Ri, for the circular solu- 
tion are marked by a filled circle, and is encircled by the Itr and 
2(T confidence ellipsoids (solid, blue lines). Corresponding values 
and confidence ellipsoids for the eccentric solution are shown with 
a triangle and dashed (red) lines. 

bust in light of the irregularities encountered with this 
data set. 



8.3. Depth and OOT Variations 

The transit depth produces a slightly excess scatter in 
reference to the global fit. We find a of 18.0 for 12 de- 
grees of freedom, which is 1.6-(T significant. The OOT is 
consistent with the global value yielding = 13.3 for 12 
degrees of freedom. As for the TDV, the varying depths 
are neither statistically significant nor particularly ro- 



9. DISCUSSION 

Due to the large number of results presented in this 
paper, we will here summarize the key findings not avail- 
able in the original discovery papers: 

■ Secondary eclipse of Kepler- 7b is detected to 3.5- 
tr confidence with depth {Fpd/Fi,) = 47 ± 14 ppm, 
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Fig. 26. — Upper Left: TTV (squares) and TDV (diamonds) for Kepler-8b (see ^2.3.1 for details). Upper Right: TTV periodogram 
(solid) and TDV periodogram (dashed) for Kepler-8b, calculated using a sequential F-test (see ^^2.3.2 for description). Lower Left: 
Transit depths from individual fits of Kepler-8b (see %2.3.1 for details). Lower Right: "Phasing" of Kepler long cadence photometry for 
Kepler-8b (see §2.5.^ for description). 



TABLE 11 

Mid-transit times, transit durations, transit depths and out-of-transit (normalized) fluxes for Kepler-8b. 
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indicative of a geometric albedo of Ag = (0.38 ± 
0.12). 

■ A marginally significant orbital eccentricity is de- 
tected for the Neptune-mass planet Kepler-4b to 
2-cr confidence, with an eccentricity of e = (0.25 ± 
0.12). 

■ A marginally significant (l.S-fi) secondary eclipse 
is detected for Kepler-5b with a depth of (26 ± 
17) ppm, for which the most plausible explanation 
is reflected light due to a geometric albedo of Ag — 
(0.15 ±0.10). 

■ A 2.6-a significance peak in the TTV periodogram 
of Kepler-6b is detected, which is not easily ex- 
plained as an alias frequency. Perturbing planets 
and exomoons are unlikely to be responsible either 
and currently our favored hypothesis is one of stel- 
lar rotation. 

■ We derive significantly different impact parameters 
for all of the Kepler planets except Kepler- 7b. 

9.1. Comparison of three hot-Neptunes 

The eccentricity of Kepler-4b is marginally detected to 
2-CT confidence. If confirmed, this would mean that three 
of the four hot-Neptunes discovered through the transit 
method have eccentricities around 0.15-0.20, the other 
two being GJ 436b and HAT-P-llb. HAT-P-26b is the 
fourth transiting hot-Neptune, recently discovered by ? 
and this system also possesses a marginal eccentricity of 
e = 0.124± 0.060. With HAT-P-26b possessing a FAP of 
12% and Kepler-4b being 11.8%, it would seem likely that 
at least one of these two systems is genuinely eccentric. 

In all cases, the tidal circularization time is expected to 
be much lower than the age of the system for Jovian-like 
tidal dissipation values, which raises the question as to 
why these planets still retain eccentric orbits. Generally 
two hypotheses have been put forward: i) an unseen per- 
turbing planet pumps the eccentricity of the hot Neptune 
ii) hot Neptunes have larger Qp values than expected. 

At the time when only one examp le of an eccen- 
tric Neptune was known, GJ 436b, Ri bas et all (|2008f ) 
made the reasonable deduction that is was more likely 
an unseen perturber was present. However, given that 
three such planets are now known, Occam's razor seems 
to favor the alternative hypothesis. For all three hot- 
Neptunes to have perturbing planets which all remain 
hidden from detection and yet produce nearly identical 
eccentricity pumping levels appears to be a less likely 
scenario than that of a common solution due to intrinsi- 
cally different tidal dissipation values and formation his- 
tory. For Kepler-4b, the observed eccentricity indicates 
Qp > 10^ if no eccentricity pumping is occurring. 

We also point out that the physical properties of 
Kepler-4b differ greatly between the eccentric and circu- 
lar fits (see Table [12]), due to the large impact of eccen- 
tricity on the lightcurve derived stellar density (jKippinel 
l2010al) . Further RV measurements are likely required to 
resolve the eccentricity of this system, since an occulta- 
tion is generally not expected to be observable given the 
very low iZp/i?*. 



9.2. Secondary eclipses 

We detect a secondary eclipse for the bloated Kepler- 
7b of depth (47 ± 14) ppm and 3.5-a confidence. This is 
above our formal detection threshold of 3-f7 thus repre- 
sents an unambiguous detection. Thermal emission is an 
unlikely source for the eclipse since the depth indicates 
a brightness temperature of 2570lg5^ which is far in 
excess of the equilibrium temperature of the planet of 
(1554 ± 32) K. A geometric albedo of Ag = (0.38 ± 0.12) 
offers a more plausible explanation for this eclipse depth. 
However, we do note that a significant nightside flux is 
apparently present which is not consistent with the re- 
flection hypothesis. We believe that this nightside flux is 
likely an artifact of the Kepler pipeline whose effects may 
mimic a long-cut filter. Although Kepler- 7b exhibits sim- 
ilarities to HD 209458b in terms of its very low density, 
the albedo clearly marks th e planet at distinc t given that 
A „ < 0.17 for H D 209458b f Rowe et al.|[20Q8l ). The study 
of lBurrows et a l. ( 2008) seems to indicate that the albedo 
requires the presence of refiective clouds, possibly com- 
posed of iron and/or silicates, as seen in L-dwarf atmo- 
spheres. It is interesting to note that Kepler- 7 is 1.3 ±0.2 
times more metal- rich than HD 209458, and the planet 
has a higher equilibrium temperature of (1565 ± 30) K 
than that of HD 209458b with (1130 ± 50) K. 

We also note that, to our knowledge and if con- 
firmed, the above measurement is the first determi- 
nation of a transiting planet's albedo at visible wave- 
lengths, and thus the first detection the associated 
reflected light. This compliments the recent polari- 
metric detection of reflected light of HD 189733b was 
made by Ber dvugina et al.l (12003) . and recently con- 
firmed in iBerdvugina et all ()2011[ ). Visible secondary 
ecli pses of other planet s have been mad e for HAT-P- 
7b (iBorucki et al. I [20091) and CoRoT-lb ([Snellen et al.l 
I2009D but in both cases the eclipse can be explained 
through a combination of reflected light plus thermal 
emission or even simply pure thermal emission for mod- 
els of large day-night contrasts, which is in fact expe cted 
for these extremely hot- Jupiters ([Fortnev et al.|[2008() . In 
contrast, thermal emission cannot explain the Kepler- 7b 
secondary eclipse unless the planet has an internal heat 
source with a luminosity equivalent to that of a M6.5 
dwarf star, which is highly improbable. 

We also flnd weak signiflcance of a secondary eclipse 
for Kepler-5b of confidence. Further transits will 

either confirm or reject this hypothesis but we note that 
the obtained depth is quite reasonable for this planet. 
Whilst thermal emission again seems unlikely based upon 
the required temperature of 2500 K versus the equi- 
librium temperature of 1800 K, a geometric albedo of 
Ag — (0.15 ± 0.10) offers a satisfactory explanation for 
the eclipse. 

For Kepler-6b and Kepler-7b, we find no evidence of 
a secondary eclipse, but the measurements do constrain 
the geometric albedos to Ag < 0.32 and Ag < 0.63 at 
3-(T confidence. For Kepler-4b, the measurement place 
no constraints on the geometric albedo. 

9.3. Differences with the discovery papers 

For almost all cases, we find significantly different im- 
pact parameters from the discovery papers, but self- 
consistent between methods A and B. There are also sev- 
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TABLE 12 

Comparison of tw o known hot-Neptunes with Kepl er-4b. GJ 43 6b values taken from VTorreA except eccentricity which comes from 

\Gillon et al.\ 1200'A) : HAT-P-llb values taken from \Bakos et al\ (20 OS) . The classification of Kepler- 4b depends greatly upon whether the 

system is confirmed as maintaining an eccentric orbit or not. 
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0.843 


0.687 
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eral examples of other parameters being different. These 
can be easily reviewed by consulting Tables [21 SI [6l [8] 
and 1101 For example, for Kepler-4b, we find a very low 
stellar density implying a larger stellar radius and thus 
larger planetary radius. Other examples of differing pa- 
rameters include the eccentricity and secondary eclipse 
parameters already discussed. 

9.4. On the error budget of planetary parameters 

The exquisite lightcurves measured by Kepler lead 
to very low uncertainties on certain parameters, even 
in the long-cadence mode. One example is the fitted 
ratio-of-radii p, where the error is about 0.5% (e.g., 
p = 0.0808 ± 0.0005 for Kepler-7b). The median error 
for the first 70 known transiting exoplanets (TEPs) on 
Rp/Ri, is 0.9%, about double that of Kepler- 7b. Limiting 
factors on the precision of this purely geometrical factor 
are the limb-darkening values for the Kepler bandpass 
and the slight degeneracy with other parameters that are 
affected by the long-cadence binning, even if re-sampled 
models are fitted (such as the b impact parameter). Us- 
ing short-cadence data, deriving precise limb-darkening 
coefficients, and accumulating many transits will signifi- 
cantly improve errors in 

The typical errors on the period for the first five Kepler 
planets are 0.00004 days, about 6 times greater than the 
median error of the period of the first 70 published TEPs, 
or abo ut 20-times the error quoted in e.g.,Hartm an et alJ 
(|2009D for HAT-P-12b, an example of a TEP discovery 
using photometry data spanning 2 years. This error on 
the Kepler planet periods, however, comes from a dataset 
with short time-span (44 days), while the ground-based 
results come from multi-year campaigns. The median 
error on the ephemeris (Tc) for the current 5 Kepler 
planets is 0.000145 days, about half that of the known 
TEPs (0.0003 days). Precision of the ephemeris data will 
greatly improve during the lifetime of the Kepler mis- 
sion just by accumulating more data, and switching to 
short-cadence mode. 

At first glance perhaps somewhat surprisingly, the er- 
rors on the Kepler planetary masses and radii are not 
different from the rest of the transiting exoplanet pop- 
ulation. The typical error on planetary masses for 70 
transiting exoplanets is 0.06 Mj, whereas the error for 
radii is 0.05 In our analysis the errors for plane- 
tary masses and radii for Kepler planets have roughly 
the same values. 

There are multiple reasons behind this. Planetary 
mass and radius scale with the respective parameters 
of the host star. In our analysis. Mi, and iZ* are de- 



termined from stellar isochrones, using Toff, [Fe/H] and 
a/Ri,. Each of these quantities have significant errors. 
The stellar T^s and [Fe/H] come from SME analysis of 
the spectra, which, in turn, due to the faintncss of the 
if ep/er targets, are of low S/N, and have larger than usual 
errors. The other input parameter, a/R^,, is related to 
the lightcurve parameters C/-R*j ^> and the RV parame- 
ters k and h. While C/^* is quite precise for the Kepler 
transits (median error 0.02 as compared to 0.15 for the 
known TEPs), the b impact parameter from the current 
LC Kepler data-set have errors comparable to ground- 
based transits. The main limiting factor, however, are 
the k and h Lagrangian orbital parameters from the RV 
data. Due to the faintness of the targets, the orbital fits 
have considerable errors in k and h. As a result, the error 
on a/i?* for the current Kepler data in this work (median 
error 0.26) is similar to that of the known TEPs (median 
error 0.2). This effect is well visible on the isochrone 
figures (Fig. [p [IHl [H [20l EH), where the 1-cr and 2-a 
confidence ellipses are shown for all Kepler planets for 
both the circular and eccentric solutions, with the latter 
clearly covering a much larger area on the isochrones. Fi- 
nally, the error in the planetary mass also scales with the 
K RV semi-amplitude, which has typical errors for the 
Kepler planets similar to the median error of the known 
planets {AK « 4ms~^). Given the faintness of the Ke- 
pler targets this is an impressive result. 

A breakthrough in both accuracy and precision of Mp 
and Rp can be expected from exploiting the extraordi- 
nary precision of Kepler, and nailin g down the stellar pa- 
ramet er s via asteroseismology (see iMulet-Marquis et alJ 
(|2009( ). iChristensen-Dalsgaard et alJ (|2010[ )). Alterna- 
tively (or simultaneously), if parallaxes for the Kepler 
targets become available (possibly from the Kepler data) , 
then anal yses sim ilar to that performed for HAT-P- 
llb ( Bak os et all [20091 can be performed, where the 
"luminosity-indicator" in the isochrone fitting will be 
the absolute magnitude of the stars, as opposed to the 
a/Ri, constraint. A final possibility is a dynamical mea- 
surement of the stellar mass by measuring the transit 
timing variations of two TEPs within the same system 
(|Agol et al.ll2005l) . 

9.5. Transit timing with LC data 

After the discovery of an exoplanet in Kepler's field, 
the subsequent photometry is switched to short-cadence 
(SC) mode. It was generally expected that this would be 
the only way to produce meaningful TTV studies. How- 
ever, we have found that by carefully correcting for the 
long integration time, the long-cadence (LC) photometry 
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yields transit times consistent with a linear ephemeris 
and of r.m.s. scatters of 10-20 seconds are possible. 

Despite the evident timing precision possible with the 
LC photometry, some issues require further investiga- 
tion. The effect of phasing, as defined in §2.3.4 remains 
somewhat unclear but in numerous instances we find that 
peaks in the TTV and TDV periodograms occur near to 
the phasing period. Subsequent aliasing and frequency 
mixing also seems to lead to numerous false positives in 
the periodograms and these issues may be diminished by 
using SC data, but are unlikely to completely disappear. 

Out of all of the studied planets, we found only one 
peak in the timing periodograms which does not ap- 
pear spurious. Kepler-6b exhibits a 2.6-cr peak of period 
(17.27 ± 0.84) days and amplitude (19.7 ± 5.0) seconds. 
We investigated the plausibility of this signal being a 
perturbing planet or exomoon but found no convincing 
supporting evidence, although such hypotheses could not 
excluded. Our favored hypothesis is that o f stellar rota- 
tion in ducing such a signal, as reported bv lAlonso et al.l 
(|2008D for CoRoT-2b. The host star has a rotational 
period of 23.5^5^9^ days and thus it is consistent with 
the TTV per i od. A bisector analysis, as performed by 
lAlonso et al.l ()2008f ). is not possible for this data due to 
the long cadence. However, the SC data will be able to 
either confirm or reject this hypothesis. 

Finally, we note that the uncertainties in all of the 
parameters from individual fits show evidence for being 
overestimates, based upon the observed scatter of these 
parameters between various transits. Both methods em- 
ployed in this analysis found very similar errors and thus 
the reason for consistent overestimation of the errors re- 
mains unclear. Both methods employed Markov Chain 
Monte Carlo (MCMC) techniques for obtaining these er- 
rors, which may be linked to the overestimation. If the 
MCMC trials are unable to sample the true global min- 
imum, perhaps due to slightly erroneous limb darkening 
coefficients or insufficient numerical resolution in the in- 
tegration time compensation procedure, the errors are 
expected to be overestimated. Nevertheless, in searching 
for evidence of signals, the accurate estimation of error 
bars is just as important as the best-fit value. We have 
proposed a F-test periodogram which is insensitive to 
absolute errors, only their relative weightings, and also 
penalizes overly-complex models. Further detailed inves- 
tigation with future Kepler timings will be possible and 
should shed light on this issue. 
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9.6. Constraints on planets, moons & Trojans 

We find no convincing evidence for perturbing planets, 
moons or Trojans in any of the systems studied here. We 
are able to exclude the presence of mean-motion resonant 
(MMR) planets for Kepler-5b, Kepler-6b and Kepler-8b 
of masses > O.79M0, OMMf^ and 0.50M® respectively, 
to 3-(T confidence. The corresponding resonant configura- 
tions are period ratios of 3:2, 4:3 and 4:3. We do not yet 
possess sufficient temporal baseline to investigate MMR 
planets for Kepler-4b and Kepler- 7b. 

Extrasolar moons at a maximum orbital separation 
are excluded by the TTV measurements for Kepler-4b 
through to Kepler-8b of masses > 11. OM®, IO.6M0, 
4.8M0, 2.5M0 and 2.IM0 respectively, to 3-(t confi- 
dence. Extrasolar moons at a minimum orbital separa- 
tion are excluded by the TDV measurements of masses 

> 13.3Me, 17.3Me, 8.2Me, 5.9Me and 21. 5M© respec- 
tively, to 3-f7 confidence. 

For Kepler-4b, Kepler- 7b and Kepler-8b, we do not yet 
possess sufficient temporal baseline to search for Trojans 
through TTV measurements. However, for Kepler-5b 
and Kepler-6b we are able to exclude Trojans of angu- 
lar displacement ~ 10° from L4/L5 of cumulative mass 

> 3.14M0 and O.67M0 respectively, to 3-a confidence. 
For all five planets, we can inspect the photometry at 
±P/6 from the transit center for signs of a photometric 
dip due the occultation of the host star by the Trojans. 
This search yielded no detections but does exclude Tro- 
jans of effective radn > 1.22i?0, 1.13i?0, O.87i?0, l.lli?0 
and 1.58i?0 for the five planets sequentially, to 3-a con- 
fidence. 
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TABLE 13 

List of important parameters used in this paper, in order of appearance. 



Parameter Name Definition 

Ml Linear limb darkening coefficient Linear term of the quadratic description of the stellar limb darkening 

U2 Quadratic limb darkening coefficient Quadratic term of the quadratic description of the stellar limb darkening 

Qc Mid companion-star separation Companion-star separation in units of stellar radii at the approximate moment of mid-1 

g Companion-star separation Companion-star separation in units of the stellar radius 

e Eccentricity Orbital eccentricity of the companion's orbit 

u) Argument of pericentre Argument of pericentre of the companion's orbit 

P Period Orbital period of the companion 

p Ratio-of-radii Ratio of the companion's radius to the stellar radius {Rp/Rt,) 

R p Radius of the companion Radius of the companion 

/?* Radius of the host star Radius of the host star 

b Impact parameter Approximate value of S when dS/dt = 0; also given by 6 = a/ R^Qc cos i 

S Sky-projected separation Sky-projected separation of the companion's centre 

and the host star's centre in units of stellar radii 

i Inclination Orbital inclination of the companion's orbit relative to the line-of-sight of the observer 

a/Rt, Semi- major axis Semi- major axis of the companion's orbit in units of the stellar radius 

a Semi-major axis Semi-major axis of the companion's orbit 

E Epoch of mid-transit Mid-transit time of the first transit in the sequence 

tc Mid-transit time Mid-transit time for a given epoch 

k Lagrange parameter k k = e cos uj 

h Lagrange parameter h k = e sin ui 

T Transit duration Time for companion to move across the stellar disc with entry 

tj First contact Instant when S = 1 + p and dS/dt < 

tji Second contact Instant when S = 1 — p and dS/dt < 

Iq Mid-transit time Instant when dS/dt = i.e. inferior conjunction 

tjij Third contact Instant when S = 1 — p and dS/dt > 

tiv Fourth contact Instant when S = 1+p and dS/dt > 

ri_4 Total duration Time for companion to move between contact points I and IV 

r2_3 Full duration Time for companion to move between contact points II and III 

Ti 2 Ingress duration Time for companion to move between contact points 1 and 2 

4 Egress duration Time for companion to move between contact points 3 and 4 

5i,4 Total eclipse duration Analogous to T\^4, but for the secondary eclipse 

T Ingress/Egress duration For circular orbits, ti2 = t34 = t 

T\ Ti duration Approximate expression for T\ from JQgging (2010a) 

p* Stellar density Average density of the host star 

C/R* Zeta over Rt The reciprocal of one half of the transit duration using the ^ingle^_&_Sackett (2005) af 

T / Rt Upsilon over /?« The reciprocal of one half of the transit duration using the j<i££ing (2010a) approximal 

Foot Primary flux Constant flux term for primary-portion of phase curve 

Foos Secondary flux Constant flux term for secondary-portion of phase curve 

Foom Middle flux Constant flux term for middle-portion of phase curve 

Fobs{t) Observed flux Flux observed by the instrument as a function of time 

_Fniodel(*) Model flux Flux from a given model as a function of time 

Ff, Stellar flux Flux received from the star 

_Fp;j Day-side flux Flux received from the day-side of the planet 

eptj Relative dayside flux Ratio of flux emitted from dayside of planet to stellar flux 

Epn Relative nightside flux Ratio of flux emitted from nightside of planet to stellar flux 

Ft Transit model flux Model flux for the primary transit 

Fg Secondary model flux Model flux for secondary eclipse 

B Blending factor Parameter to quantify amount of blended light in the aperture 

<I> Phase function Flux received from planet+star-|-blend without transits or eclipses 

7rei RV offset Constant used for offsetting the radial velocity measurements 

K RV semi-amplitude Semi-amplitude of the radial velocity signal 

T^ff Effective temperature Effective temperature of the star 

[Fe/H] Metallicity Metallicity of the host star 

V sin i Projected rotation velocity Projected stellar rotation velocity 

Mt Stellar mass Mass of the host star 

Mp Planetary mass Mass of the planet 

L» Luminosity Stellar luminosity 

AIv Absolute magnitude Absolute magnitude of the host star 

pp planetary density Density of the planet 

g Stellar surface gravity Surface gravity of the host star 

gp Planetary surface gravity Surface gravity of the planet 

Teq Equilibrium temperature Equilibrium temperature of the planet 

Tcirc Circularization timescale Time for eccentricity to reduce by an e-fold 

Qp Quality factor of planet Quality factor of planet due to tidal dissipation 
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TABLE 14 

Comparison oj various models tried for the radial velocity data of each system, in method A. Results come from a preliminary MCMC 
plus simplex. Details regarding the differences between the various models is provided in %2.1. The Bayesian Information Criterion 
(BIC) is used to select the preferred model for the circular and eccentric fits, which is then used in the final analyses. 



Model 


Circular Orbit BICs 


Eccentric Orbit BICs 


Kepler-4 


Free 7; Free Aitroj 


34.395 


36.354 


Free 7; Fixed Atnoj 


31.144 


33.300 


Fixed 7; Free Attroj 


31.486 


35.732 


Fixed 7; Fixed Attroj 


29.556 


32.685 


Kepler- 5 


Free 7; Free Attroj 


29.275 


29.222 


Free 7; Fixed Attroj 


26.768 


28.487 


Fixed 7; Free Attroj 


43.990 


47.592 


Fixed 7; Fixed Attroj 


45.321 


45.889 


Kepler- 6 


Free 7; Free Attroj 


17.467 


19.790 


Free 7; Fixed Attroj 


21.234 


17.693 


Fixed 7; Free Attroj 


15.592 


17.878 


Fixed 7; Fixed Attroj 


20.157 


16.834 


Kepler- 7 


Free 7; Free Attroj 


13.195 


17.276 


Free 7; Fixed Attroj 


13.0113 


15.518 


Fixed 7; Free Attroj 


11.508 


15.530 


Fixed 7; Fixed Attroj 


11.160 


13.664 


Kepler-8 


Free 7; Free Attroj 


30.877 


37.213 


Free 7; Fixed Attroj 


28.162 


34.311 


Fixed 7; Free Attroj 


28.163 


34.340 


Fixed 7; Fixed Attroj 


27.130 


31.675 



TABLE 15 

Outliers identified during our analysis for each system. 



(BJD-2454900) 


Standard deviations 


(BJD-2454900) 


Standard deviations 


Kepler-4 


7.S.176 


3.5 


93.590 


5.3 


81. ()(>() 


-5.1 


90.1)82 


- *") 7 


Kepler- 5 


55.418 


3.9 


53.641 


4.1 


62.795 


5.3 


54.437 


9.6 


70.683 


4.3 


71.030 


-3.8 


77.753 


3.7 






Kepler- 6 


61.733 


4.1 


56.236 


4.1 


67.945 


-3.5 


66.024 


3.6 


86.070 


3.9 


88.011 


4.2 


55.929 


3.5 


91.775 


1.0 


Keptcf-V 


72.113 


-3.6 


84.415 


-3.5 


96.287 


-3.9 


94.571 


3.6 


Kepler-8 


78.192 


10.1 


56.656 


3.7 


58.883 


17.3 


58.781 


4.2 


60.681 


4.1 


73.349 


3.7 


.I.'IO:-!!) 


1.7 







